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^ ! ABSTRACT 

' We report on the analysis of a suite of SPH simulations (incorporating eooling and 

, star formation) of mergers involving idealised X-ray clusters whose initial conditions 

' resemble relaxed clusters with cool compact cores observed by Chandra and XMM. 

00 ' The simulations sample the most interesting, theoretically plausible, range of impact 

, parameters and progenitor mass ratios. We find that all mergers evolve via a common 

progression. We illustrate this progression in the projected gas density, X-ray surface 
I brightness, SZ, temperature, and gas entropy maps. Several different classes of tran- 

r~| . sient "cold front"-like features can arise over the course of a merger. Each class is 

i' distinguished by a distinct morphological signature and physical cause. We find that 

' I all of these classes are present in Chandra and XMM observations of merging systems 

^ . and propose a naming scheme for these features: "comet-like" tails, bridges, plumes, 

' streams and edges. In none of the cases considered do the initial cool compact cores 

, of the primary and the secondary get destroyed during the course of the mergers. In- 

1^' ' stead, the two remnant cores eventually combine to form a new core that, depending 

on the final mass of the remnant, can have a greater cooling efficiency than either of 
its progenitors. We quantify the evolving morphology of our mergers using centroid 
;J] ■ variance, power ratios and offset between the X-ray and the projected mass maps. We 

I find that the centroid variance best captures the dynamical state of the cluster. It also 

provides an excellent indicator of how far the system is from virial and hydrostatic 
equilibrium. Placing the system at z = 0.1, we find that all easily identified observable 
traces of the secondary disappear from a simulated 50 ks Chandra image following the 
second pericentric passage. The system, however, takes an additional ~ 2 Gyr to relax 
and virialize. Observationally, the only reliable indicator of a system in this state is 
the smoothness of its the X-ray surface brightness isophotes, not temperature fiuc- 
tuations. Temperature fiuctuations at the level of AT/T 20%, can persist in the 
final systems well past the point of virialization, suggesting that that the existence of 
temperature fluctuations, in and of themselves, do not necessarily indicate a disturbed 
or unrelaxed system. 

Key words: cosmology: theory - galaxies: clusters: general - intergalactic medium 
- X-rays: general 



1 INTRODUCTION 

The utility of galaxy clusters as tools for studying cosmol- 
ogy and the history of structure formation is well estab- 
lished. Furthermore, sizable samples of clusters extending 

* E-mail: gbpoole@uvic.ca 



to high redshift will soon be available from wide field X- 
ray and Sunyaev-Zel'dovich (SZ) surveys, greatly expanding 
their importance and ensuring that they will feature promi- 
nently in these pursuits throughout the foreseeable future. 

Often implicit in such studies is the assumption that 
clusters are in equilibrium. However, galaxy clusters are 
dynamically young systems and are thus prone to regu- 
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lar disruptions by mergers JCohn fc White! l2005t) . Recent 
Chandra and XMM observations support this. High res- 
olution spectroscopic imaging has provided dramatic ev- 
idence of ongoing mergers in several systems. As a re- 
sult, we now know that the intracluster medium (ICM) 
of merging clusters can express a varie ty of complicated 
trans ient features i ncluding cold fronts ijM arkevit ch et alJ 
l200d) . shock fronts jMarkevitch fc Vikhlininll200ir) and the 
"sloshing" of cool cores in their dark matter potentials 
jMarkevitch. Vikhhnin. fc Mazzottal I2001I : iDupke fc White! 
!2003^ . 

Such observations have stimulated several theoretical 
and computational studies de voted to und e rstan ding the 
impact of mergers on clusters. !Bialek et alj l!2002ll , for in- 
stance, have shown that cold fronts can be produ ced throug h 
the di sruption of cold cores during mergers while lMotl et alJ 
||200^ have shown that a significant fraction of a merging 
cold core can avoid shock heating, contributin g to the assem- 
bly of the remnant's cold core. Furthermorer iRowlev et'all 
i2004l) have shown that smooth accretion and mergers lead 
to different evolutions in a system's X-ray l uminosity and 
tempe rature. Combined with the results of !0'Hara et alJ 
(I2OO3) who suggest that all clusters may exhibit departures 
from equilibrium due to the lingering effects of their past 
merger activity, it becomes clear that we must understand 
the impact of merging in accounting for the observed prop- 
erties of clusters; even the seemingly relaxed ones. 

Determining the detailed effects of a single merger on 
a cluster evolving in a cosmological environment can be 
greatly complicated by the presence of extant substructure, 
subsequent mergers and smooth accretion. Simulations of 
idealised two-body cluster mergers provide an alternative 
approach and, given the fact that most clust ers obtain much 
of th eir mass through significant mergers JCohn fc White! 
[200^, are not a significant departure from cosmological sim- 
ulations. They allow control of the initial structure and or- 
bits of the interacting systems and facilitate the investiga- 
tion of not only the obvious short-term effects of the event, 
but more subtle long-term effects as well. They also simplify 
the interpretation of specific merger observations which typ- 
ically involve only two dominant components. 

Seve ral authors have explored cluster mergers in this 
wav fsee !Rqe ttiger eraL 1993; Schindlcr fc MucUei! [l993l: 
iPearce etaLn994 : iRoettiger e ^,L.199.7., for pione ering stud- 
ies). Amoung recent studies. ICom^^t alJ (j2003) have con- 
structed a suite of idealised merging systems of 4:1 and 16:1 
mass ratios to study the effects of mergers on the stability 
of central cool cores. The initial conditions for their clus- 
ters were chosen to have gas core radii of ~ 250 kpc and 
the dark matter in the ir simulations was distrib uted accord- 
ing to a King model. iRicker fc SarazinI ll200ll . RSOl here- 
after) have studied mergers between idealised systems with 
cuspy dark matter profiles. Since the focus of their study 
was on the luminosity and temperature "boosts" which a 
merger system experiences upon impact of the cluster cores, 
they did not include the effects of cooling in their simula- 
tions. Lastly, Ritchie fc Thomas ( 2002) have analysed a set 
of simulations incorporating cooling and star formation for 
idealised merging systems constructed with cuspy dark mat- 
ter profiles. Their initial conditions primarily examine sys- 
tems with long cooling times and core radii of Vc ~ 100 kpc, 



although three cases represent systems with short cooling 
times and small core radii (vc ~ 50 kpc). 

However, it is now known that clusters have cuspy 
dark matter profiles and that 70 — 90% of observed sys- 
tems possess com pact cores (rc ~ 50 kpc. IPeres et a'l]!l998l : 
lEdee et al.lll992^ with cooling times short enough to be of 
significant importance to their structure. Hence, a compre- 
hensive study of cluster mergers utilising initial conditions 
faithful to our contemporary understanding of a typical clus- 
ter's structure while including the important effects of cool- 
ing has not been performed. Given the wide ranging appli- 
cations of cluster studies, the importance of mergers to un- 
derstanding the structure of the ICM and the wealth of de- 
tailed cluster merger observations requiring interpretation, a 
thorough revisit of the issues introduced by cluster mergers 
would be of great utility. 

In this paper and the series to follow, we will present 
an analysis of a suite of SPH simulations of isolated two- 
body mergers between idealised clusters constructed to pos- 
sess realistic compact cool cores. With this study we seek to 
provide a conceptual foundation for subsequent studies to 
be conducted in a full cosmological context. Our approach 
utilises insights provided by recent high resolution observa- 
tions of clusters to motivate the initial conditions for a set of 
nine hydrodynamic simulations (incorporating the effects of 
cooling and star formation) involving a representative range 
of mass ratios and impact parameters. In this article we 
explain our approach to constructing idealised merger sim- 
ulations and present their general dynamical progression, 
eventual relaxation and describe the transient phenomena 
which manifest during the event. Since many of the phe- 
nomena we study have low surface brightnesses, we utilise 
synthetic observations to delineate realistically observable 
phenomena from those which are not. In subsequent papers 
we will examine in detail several specific issues which we can 
only touch upon in the present work. 

In Section |21 we describe our methods for initialising 
merging systems and placing them on realistic orbits as well 
as the details of our numerical methods. In Section |21 we 
qualitatively describe the evolution of our cluster mergers 
through a generic sequence of states and introduce the tran- 
sient structures formed in the process. In Section |1| we ex- 
amine the evolution of several measures of our mergers' ap- 
parent degree of disruption to determine when our systems 
would appear undisturbed under reasonable observational 
circumstances. In Section |S] we assess the degree to which 
our apparently relaxed merger remnants are formally so. In 
Section |^ we examine in further detail the different tran- 
sient structures which form during a merger, describing their 
properties and the processes which drive their creation. Fi- 
nally we summarise our study in Section |7| 

In all cases our assumed cosmology will be 
(!:!Af,f^A)=(0.3,0.7) with Ho=75 km/s/Mpc and 

nt = 0.02/1"^ 



2 SIMULATIONS 

We run our simu l ations with GASOLINE 
iWadslev. Stadel. fc Quimil l2004l) . a versatile parallel 
SPH tree code with multi-stepping. We include the effects 
of radiative cooling, star formation and minimal feedback 
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Table 1. Initial separations (r), radial and transverse velocities (vr and vt respectively) for our simulation orbits and 
measures of the final remnant's global structure (the spin parameter A200 and triaxial ratios (J200 and S20O1 all measured 
within -R20o)- At the start of our simulations, t = 0; when the secondary's centre of mass passes the virial radii Hioo or 
^200 of th6 primary, t = t'^ or t = to respectively. 



Mp-.Ms 




to[Gyrs] 


r{0)[Mpc] 


vt{0)[km/s] 


t)r(0)[fcm/s] 


^min [kpc] 


A200 


9200 


■5200 


1:1 


0.00 


3.4 


6.54 





-721 


5 


0.0 


0.67 


0.62 


1:1 


0.15 


3.4 


6.54 


-111 


-708 


110 


0.021 


0.70 


0.65 


1:1 


0.40 


3.2 


6.54 


-315 


-925 


390 


0.059 


0.68 


0.65 


3:1 


0.00 


3.4 


5.53 





-574 


12 


0.0 


0.76 


0.73 


3:1 


0.15 


3.4 


5.53 


-87 


-587 


103 


0.014 


0.80 


0.75 


3:1 


0.40 


2.6 


5.53 


-271 


-908 


360 


0.042 


0.78 


0.76 


10:1 


0.00 


2.2 


4.81 





-905 


10 


0.0 


0.96 


0.94 


10:1 


0.15 


2.4 


4.81 


-106 


-872 


130 


0.007 


0.96 


0.92 


10:1 


0.40 


2.2 


4.81 


-284 


-950 


419 


0.018 


0.97 


0.93 



from supernovae in our simulations but the effects of 
feedback from active galactic nuclei (AGN) are omitted. 
It is important to note that we do not include the effects 
of magnetic fields, pressure generated from cosmic rays 
or conduction and that our choice of algorithm is not 
well suited to modelling the effects of turbulence. The 
relevance of each of these processes is not well established 
(theoretically or observationally) and we can only provide 
the caveat that some of our results are susceptible to 
change once we ha ve the resources to properly model 
them (although, see iMarkevitch et al] l2003l : iFabian et alJ 
l2003l . for references suggesting that neither conduction nor 
turbulence are relevant in cluster cores). 

In this section we describe how our clusters are ini- 
tialised, how we set up pairs of clusters to form a merging 
system, the numerical methods and code parameters of our 
runs. 



2.1 Cluster initial conditions 

The structure of a merger remnant is likely to depend sen- 
sitively on the init ial structure of the interacting systems 
jGomez et al. ''2002'l making it important to implement re- 
alistic initial conditions for the gas and dark matter prop- 
erties of our cluste rs. W e follow the analytic pr escription of 
iBabul et all (|20o3) and lMcCarthv et all ^l2004^ to produce 
systems which conform with recent theoretical and observa- 
tional insights into cluster structure. 

The dark matter den sity profiles of our systems fol- 
low an NFy-hk e form jNavarro. Frenk. fc Whlt^ Il996l: 
iMoore et ai]ll998l) given by 



{r/rs)P (1 + [r/rs)) 

with the central asymptotic logarithmic slope chosen to 
be /3 = 1.4 and selected to yield a concentration c — 
R20o/rs ~ 2.6 (we will use i?A throughout to indicate the 
radius within which the mean density of the system is A 
times the critical density, pc = 3Hq /SttG). This dark mat- 
ter distribution would yield c = 4.5 if fit by Eqn. with 



/3 = 1 and is con sistent with resu lts from cosmological N- 
body simulations jEke et alJl2'00lh . 

The initial density and temperature profiles of the clus- 
ters are set by requiring that (1) the gas be in hydrostatic 
potential within the halo, (2) the ratio of gas mass to dark 
matter mass within the virial radius be flb/{0.m — ilt), and 
(3) the initial gas entropy® scale as S{r) oc r^'^ over the bulk 
of the cluster body. 

Our choice for the initial gas distribution is motivated 
by the fact that we are interested in studying "typical" sys- 
tems which correspond to the most likely pre-interaction 
configuration for merging clusters. Estimates in the liter- 
ature suggest that the majority of observed clusters (70- 
90%) possess dense, cool, compact cores with short cooling 
times (Peres et al. 1998). Recent Chandra observations re- 
veal that these systems exhibit power-law entropy profiles 
of the form S{r) oc t-I O-I-^^ beyond the central 10 kpc, over 
a wide range of mass scales jPonahue et al.ll2005ft . We note 
that our choice for the form of the initial density profile 
conforms not only with these observations but also with the 
resul t s of high resolu tion cluster simulations (|Lewis et alJ 
l2000l:IVoit et alJl2003t) . 

We normalize the entropy profiles such that the tem- 
perature of the ICM at R^ir is half the virial tempera- 
ture. Thi s also approximately matches cosmo logical sim- 
ulations llLewis et alJl200d: iLoken et alJ l2002l) and obser- 
vations (^D^^ranx^^^Motold^l]2002^ . Finallv. for practical 
purposes, we start our systems with small, constant, low- 
entropy (lOkeVcm'^) cores. The cores, however, evolve to 
equilibrium power-law distributions quickly (~ 0.5 Gyr) and 
well before the clusters begin interacting significantly. 

Particle realizations are initialised using the ZENO 
package of J. Barnes. Gas particles are distributed 
spherically and assigned temperatures according to the 
model described above. The dark matter velocity dis- 
tribution is taken to be isotropic and the distribu- 
tion function of particle energies calculated by solv- 
ing the Abel integral equation (Binnev fc Tremain3 Il98^ : 



t For colour plots, please see the electronic version or 
http: / /astro. phys.uvic.ca/~babul/Merger_PaperI 
® We use the standard proxy for entropy given by S = kT jr^J^ 
with rie and T representing the electron density and temperature 
of the gas. 
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Figure 1. Maps of the gas surface density for our 1:1 simulations, projected along an axis normal to the plane of the orbit (3Mpc 
on a side). White contours are X-ray surface brightness isophotes from simulated 50 ks z = 0.1 Chandra observations (3x to 48 X 
the background stepping by factors of 2). Times depicted are iciosest ~200 Myr (where tc.iosast'^^ the time of first pericentric passage), 
tcioaest+lOO Myr, tapo(where tapois the time of first apocentric passage), taccrete+100 Myr (where taccrete\& the time of the second 
and usually the final pericentric passage) and t^^i^x , the time when system appears relaxed. All times are measured relative to to (the 
time when the secondary's centre of mass first crosses -R200 of the primary) 



iKazantzidis. Magorrian. fc Moor3 |2004|) . Previous cluster 
merger studies have been initialised with Gaussian distri- 
butions resulting in systems which subsequently evolve sig- 
nificantly from their initial conditions (RSOl). Our approach 
initialises the system very close to equilibrium, ensuring that 
our interacting systems accurately maintain our desired ini- 
tial conditions until they collide. 

Negative values for the dark matter energy distribu- 
tion function are avoided by truncating our clusters (dark 
matter and gas) smoothly beyond the virial radius {Rvir) 
with the function p(r) oc [r / Ryir)'^ ey.]i{~[{r — Rmr) / Rsm]^) , 
where 7 is set to ensure a continuous first derivative for 
p(r). Beyond the point where the gas pressure reaches that 
of the intergalactic medium (IGM) (we choose Pigm = 
{QbPckTiGM)/{fimH) with TiGM = 3 x 10^ K), we sur- 
round our systems with a dynamically negligible uniform 
gaseous medium of pressure Pigm- This is done mainly to 
confine small numbers of high-velocity gas particles occa- 
sionally produced during our simulations. Without this ex- 
ternal medium, these particles would adiabatically expand 
to unrealistically low densities and temperatures and pose 
numerical difficulties for the code. 

In the present study we examine mergers between sys- 
tems set to have mass ratios of 1:1, 3:1 and 10:1. In what 
follows we shall distinguish the two interacting clusters by 
identifying the most massive 10^' system as the primary 
system (with mass Mp) and the smaller incident system 
(with mass Ms) as the secondary. For our equal mass merg- 



ers, we arbitrarily choose one system as the primary with no 
eflect on our results due to the symmetry in such cases. All 
positions and velocities are measured relative to the centre 
of mass of the primary system. 

The virial mass of primary system is set to 10^^ M© in 
all cases and the secondary systems have virial masses of 
10^^ Mq, 10" '^ Mq, and lO" Mq. These values are larger 
than M200, which takes the values 7.5 x 10^'' Mq, 2.4 x 
10" Mq, and 7.8 x 10^^ Mq for the three systems. We shall, 
however, refer to the systems by their virial mass rather than 
M200- The circular velocities of the halos at J?200 (V200) are 
1340 km s~^, 920kms~^, and 630kms~^ respectively. 



2.2 Initial kinematics 

The problem of initialising the kinematics of idealised 
mergi ng clusters has been faced by several authors 
(e.g. iRitchie fc Thomas! I200I RSQl). Whereas previous 
studies have initialised orbits based on analytic arguments, 
our initial conditions are motivated by analysis of the or- 
bital properties of substructure in high-resolution cosmolog- 
ical dark matter simulations. 

We construct orbits for these systems which produce 
specified radial and tangential velocities for the secondary 
system {vr and vt respectively) when its centre of mass 
reaches the virial radius of the primary (Rvir)- For each 
of the three mass ratios we study, we examine three orbits 
selected to produce a typical value of Vr(Rvir) and to cover 
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Figure 2. Maps of the projected "spectroscopic-like" temperature iMazzotta et al]|200'J) maps for a 0.5 Mpc thick shcc (3 Mpc on a 
side) passing through the centre of our 1 : 1 simulations. Contours depict the SZ effect (for the entire simulation volume along the line 
of sight) and represent log(y) = —5.5 to —4 in increments of 0.25^. The times represented are the same as Fig. ITl 




Figure 3. Entropy maps for a 0.5 Mpc thick slice (3 Mpc on a side) through the centres of our 1 : 1 simulations^ . The times represented 
are the same as Fig. 
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Figure 4. Maps of the gas surface density for our 3:1 simulations, projected along an axis normal to the plane of the orbit (3Mpc 
on a side). White contours are X-ray surface brightness isophotes from simulated 50ks 2; = 0.1 Chandra observations (Sx to 48 X the 
background stepping by factors of 2). Times depicted are tcjosest ~200 Myr, fcjoses^-l-lOO Myr, tapo, taccrete+100 Myr and tj-^i^x- See 
Figure 1 caption or Table 2 for a brief definition of the timescales. All times are measured relative to the to when the secondary's centre 
of mass crosses the virial radius of the primary for the first time^ . 



a significant range of tlie transverse velocity Vt{Rvir) giving 
rise to mergers found in cosmological dark matter simula- 
tions. 

In an N-b ody study of haloes merging onto rich clusters, 
iTormenI il997ft find that the average velocity of an accreting 
secondary system is v{Rvir) = (1.1 ± 0.1)Vc(-Ri,ir), where Vc 
is the primary s ystem 's circular velocity. In a nother study, 
IVitvitska et"all ll2002h (see also iBensonllioOsI) examine the 
distribution of secondary subhalo velocities moving within 
the virial radius of massive haloes. They find that the total 
velocities of accreting systems {i.e. those with negative ra- 
dial velocities in the range r = 0.8-1.2_R„i, ) relative to the 
primary are distributed normally. They report mean infall 
velocities of (v) ~ Vc (independent of mass), consistent with 
[Tormcn ( 1997). They also find that tangential 2-dimensional 
rms velocities decrease with the secondary mass, ranging 
from (j_L ~ 0.4Vc for 1:1 mergers to a± « 0.67Vc for 10:1 
mergers. It is important to note that these statistics likely in- 
clude systems making secondary or tertiary encounters with 
their primary. The orbits of such systems have likely been 
rendered more isotropic, making these tangential dispersions 
upper limits for substructure accreting for the first time. 
To produce systems which will merge fully within a rea- 
sonable amount of time, we have chosen to examine three 
values of vt/Vc which roughly span the lower half of the 
? distribution. Specifically, we construct orbits to produce 
Vr{Rvir) = 1.2Vc{Rvir) (slightly higher but consistent with 
the mean velocities cited above) and vt/Vc ~ 0, 0.15, and 
0.4. Unlike in the rest of our study, we define the virial ra- 



dius here by an overdensity A — 100, corresponding to the 
threshold used by ?. Doing this for mass ratios of 1:1, 3:1 and 
10:1 produces nine simulations which constitute the basis of 
the analysis we present in this study. 

Throughout our analysis we shall refer to three coordi- 
nate axes; x and y will denote directions in the plane of the 
initial orbit and z the direction orthogonal to this plane. To 
properly capture the initial effects of the encounter as the 
outermost envelopes of the systems make contact and the 
interaction shocks begin to form, we initialise our clusters 
separated along the x axis so that their smoothly truncated 
envelopes barely touch {i.e. separated by slightly more than 
the sum of their virial radii). As a result, the initial con- 
ditions we have outlined above must be imposed at a time 
part way through the run when the centre of mass of the 
secondary system crosses the virial radius of the primary. 
To determine the initial velocities which yield these condi- 
tions, we first approximate the systems as point masses and 
analytically determine the initial conditions which yield the 
desired constraints. To account for the extended, tidally de- 
forming mass distributions, we then iteratively refine these 
orbits with low-resolution dark-matter-only simulations. 

Table lists the parameters relevant to the initial con- 
ditions of our runs including the centre-of-mass separation 
of the secondary from the primary (r), initial tangential and 
radial centre-of-mass velocities of the secondary relative to 
the primary {vt and Vr), the time of virial crossing {to, using 
the virial radius -R200), and the centre-of-mass distances of 
closest approach {rmin). As a test of the validity of our or- 
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Figure 5. Maps of the projected "spectroscopic-like" temperature iMazzotta et al]|200'J) maps for a 0.5 Mpc thick shcc (3 Mpc on a 
side) passing through the centre of our 3 : 1 simulations. Contours depict the SZ effect (for the entire simulation volume along the line 
of sight) and represent log(y) = —5.5 to —4 in increments of 0.25^. The times represented are the same as Fig. Bl 




Figure 6. Maps for a 0.5 Mpc thick slice (3 Mpc on a side) through the centres of our 3 : 1 simulations^. The times represented are 
the same as Fig. El 
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Figure 7. Maps of tlie gas surface density for our 10:1 simulations, projected along an axis normal to the plane of the orbit (3Mpc 
on a side). White contours are X-ray surface brightness isophotes from simulated 50ks 2 = 0.1 Chandra observations {3x to 48 X the 
background stepping by factors of 2). Times depicted are tcjosest— 200 Myr, fcjosest+lOO Myr, tapo, taccrete+100 Myr and t^eiax- See 
Figure 1 caption or Table 2 for a brief definition of the timescales. All times are measured relative to the to when the secondary's centre 
of mass crosses the i?200 radius of the primary for the first timet 



bits we have compiled the spin parameter and dark matter 
triaxial shapes for our final combined systems. We compute 
the spin parameter from A = J\E\^''^ /GM^^^ with J being 
the total angular momentum, E the total energy and M the 
total mass of the remnant, all integrated within ii2oo- We 
compute the triaxial axes rati os at -R200 (denoted ^ 2 00 an d 
S200) following the approach of lDubinski fc Carlberd il99ll) . 
All are in reasonable agreement with accepted distributions 
from cosmological simulations l|Bajnes fc Efstathioulll98'^ 
IJing fc Sutoll2002h . 



2.3 Numerical methods 

To compute the hydrodynamic forces between gas par- 
ticles, ^^e__use__GA^OLIN^iiiits default configura- 
tion iWadslev. Stadel. fc Quiniil l2004ll : we use the Benz 
arithmetic-asymmetric implementation of the SPH momen- 
tum and energy equations. The gas quantities are smoothed 
over 32 particles, the Courant time step parameter is set 
to ?7c = 0.4, and the viscosity uses the standard Monaghan 
formulation of viscosity with a = 1 and /3 = 2. 

The gravitational forces are calculated with a tree that 
uses an opening angle — 0.8. Each particle trajectory is 
integrated on its own time step Ati, which is related to the 
particle's acceleration at by the criterion Ati < ri^e/aiY^'^ 
with tolerance parameter 77 = 0.2. The gravitational soften- 
ing uses spline kernel interpolation, and the spline softening 
length is set to a constant value of e = 10 kpc for all cluster 



particles. This ensures that gravitational forces in the core 
are well resolved. 

Radiative cooling allows for the formation of resilient 
dense gaseous cores. In some cases it can also dissipate sig- 
nificant amounts of the thermal energy deposited in the cen- 
tres of merger remnants, further stabilising them. Since we 
are acutely interested in the evolution of these regions for 
this series of papers, the inclusion of cooling in our simula- 
tions is necessary. 

We focus our analysis on radii r > 40 kpc ~ 0.03i?5oo 
and have sought to ensure that we accurately model the 
behaviour of the system in that range. This is because the 
precise nature of very central regions of the observed com- 
pact cool cores is still not well unde rstood. The tempera - 
ture floors observed in their centres JPeterson et alJl2003h 
indicate the presence of a source of heating whose nature 
is still highly uncertain with several candidate mechanisms 
presently being considered (e.g. heating from AGN, con- 
duction and turbulent mixing to name only a few). There 
is also evidence of other complex phenomena occurring in 
the central 40 k pc, such as multiphase cooling and ongoing 
star formation jjaffe et al . 2005: Ega mi et al . 200fi) . These 
processes cannot be realistically captured with SPH simula- 
tions such as ours. In light of these issues, we shall compute 
"central" quantities at 40 kpc and shall exclude the central 
30 kpc ~ 0.025-R500 of the system from integrated quantities 
such as luminosities and globally averaged temperatures. 

For our simulations we use a cooling function appropri- 
ate to a primordial gas. However, typical clusters have metal- 
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Figure 8. Maps of the projected "spectroscopic-like" temperature ("Mazzotta et al."'2004') maps for a 0.5 Mpc thick shcc (3 Mpc on a 
side) passing through the centre of our 10 : 1 simulations. Contours depict the SZ effect (for the entire simulation volume along the 
line of sight) and represent log{y) = —5.5 to —4 in increments of 0.25^. The times represented are the same as Fig. 171 
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Figure 9. Maps for a 0.5 Mpc thick slice (3 Mpc on a side) through the centres of our 10 : 1 simulations^. The times represented are 
the same as Fig. 171 
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Table 2. Times of important stages in the dynamical evolution of our simulations measured relative to to, the time when the 
secondary's centre of mass crosses the primary's -R200 radius. These are as follows: tdosest is the time of first pericentric passage; 
tapo is the time of the first apocentric passage; taccrete is the time of the second pericentric passage typically during which the 
secondary's gas core merges with that of the primary and it also marks the beginning of stream accretion; t^f,i^^ is the time at 
which the total system visually appears relaxed in simulated 50 ks Chandra observations at z = 0.1; tyi^i^l is time at which the 
system's dark matter is assessed to be virialised within -R500; ^■nd thydro is the time at which the system returns to hydrostatic 
equilibrium (to within 10%) at -R500 . Quantities in brackets indicate the redshift at which the secondary core would have to cross 
R200 of the primary for the system to be observed in the given state at 2 = 0.1 and z = 0.5 respectively. Additional parameters in 
the table are: the mass ratio of the primary to the secondary (column 1) and the secondary's transverse velocity at when its centre 
of mass crosses the virial radius of the primary, in units of the primary's circular velocity (column 2). 



Mp-.Ms 


Vt/Vc 




^closest 




^accrete 




^virial 


^hydro 


1:1 


0.00 


0.5 


(0.14, 0.58) 


1.2 (0.21, 0.70) 


2.2 (0.33, 0.92) 


4.4 (0.65, 1.73) 


4.4 (0.65,1.73) 


7.9 (1.73,>10) 


1:1 


0.15 


0.6 


(0.15, 0.60) 


1.5 (0.25, 0.77) 


2.3 (0.34, 0.95) 


4.4 (0.65, 1.73) 


4.4 (0.65,1.73) 




1:1 


0.40 


0.5 


(0.14, 0.58) 


1.4 (0.24, 0.75) 


2.7 (0.39, 1.06) 


5.4 (0.85, 2.44) 


3.8 (0.55,1.44) 


3.3 (0.47,1.25) 


3:1 


0.00 


0.5 


(0.14, 0.58) 


1.4 (0.24, 0.75) 


2.3 (0.34, 0.95) 


5.0 (0.77, 2.10) 


3.4 (0.49,1.28) 


4.3 (0.63,1.67) 


3:1 


0.15 


0.5 


(0.14, 0.58) 


1.5 (0.25, 0.77) 


2.5 (0.26, 1.00) 


4.6 (0.69, 1.84) 


3.4 (0.49,1.28) 


4.3 (0.63,1.67) 


3:1 


0.40 


0.8 


(0.18, 0.63) 


1.9 (0.29, 0.85) 


3.5 (0.50, 1.32) 


5.7 (0.92, 2.75) 


4.3 (0.63,1.67) 


3.6 (0.52,1.36) 


10:1 


0.00 


0.8 


(0.18, 0.63) 


2.2 (0.33, 0.92) 


3.5 (0.50, 1.32) 


4.8 (0.73, 1.96) 


4.6 (0.69,1.84) 


5.2 (0.81,2.26) 


10:1 


0.15 


0.6 


(0.15, 0.60) 


2.1 (0.32, 0.90) 


3.8 (0.55, 1.44) 


5.6 (0.90, 2.65) 


4.5 (0.67,1.78) 


5.1 (0.79,2.18) 


10:1 


0.40 


0.8 


(0.18, 0.63) 


2.5 (0.26, 1.00) 


4.8 (0.73, 1.96) 


7.8 (1.68,>10) 


5.6 (0.90,2.65) 


5.8 (0.95,2.88) 



licities of Z ~ O.SZq and by omitting heavy element line 
emission, we are underestimating the cooling rate of mate- 
rial with temperatures T — 10^-10^ K. A very small fraction 
(< 1%) of the gas mass in our systems exists at these temper- 
atures; most of the mass has T = lO'^-lO^ K. Additionally, 
this discrepancy has a negligible effect on the dynamics of 
our simulations. All subsequent analyses {e.g. producing X- 
ray surface brightness maps) are comput ed assuming a more 
realis tic Raymond and Smith model iRavmond fc SmithI 
Il977l) compiled for a metallicity of Z ~ O.SZq. 

We include star formation in our simulations using an 
algorithm which turns cold and dense gas into collisionless 
star p articl es, using a recipe similar to that of iKatz et alJ 
il996l) and IStinson et a l. (2006). We first select gas parti- 
cles which are capable of forming stars on the basis of four 
eligibility criteria: the gas must be dense {nn > 0.1 cm~'^); it 
must be either cool (T < 3 x lO'* K) or cooling {dT/dt < 0); 
its flow must be converging (V-w < 0); and the particle must 
be Jeans unstable. Once selected, star formation produces 
star particles with a probability p = 1.0 — exp(— c*^r/f/orm) 
where the star-formation efBciency is taken to be c, — 0.1. 
We set tform to the maximum of the cooling time and the 
dynamical time, unless the particle is already cool (T < 
3 X 10* K) in which case we always set it to the dynami- 
cal time. If the particle is chosen to form stars, a new star 
particle is created with 1/3 of the mass of a full-sized gas 
particle. 

Star formation and supernova explosions go hand in 
hand. The energy released by the supernovae is injected into 
the surrounding intracluster medium as follows. At every 
timestep, we consider each star particle in the simulation 
and based on its age, determine the associated supernovae 
rate. The total energy released by the supernovae over the 
course of the timestep is then injected at a steady rate to the 
surrounding gas particles over the timestep. The star forma- 
ti on and feedback reci pe we use is similar to that described 
in lStinson et al.l i2006l) . except that we do not disable radia- 
tive cooling during feedback, and stellar winds which release 
mass from intermediate-mass stars are not included. This is 



a "minimal" feedback scheme, since the energy is quickly ra- 
diated away and the feedback does not impact the evolution 
of the ICM in any meaningful fashion. 

Finally, we note that our chosen mass resolution models 
the initial conditions of a 10^^ Mq system with 217441 dark 
matter and 211914 gas particles, where about 70% of the 
particles in each species are within i?2oo- The dark matter 
particle mass is 4.4 x 10^ M0 and the initial gas particle mass 
is 5.9 x 10* M0. In test runs, we have found that the cooling 
rates of isolated clusters are not robust until there are ~ 100 
particles within the "cooling radius" (i.e. the radius within 
which tcooi < t Hubble; Rc ~ 150 kpc for our simulations). 
With our initial conditions, the central core (rc ~ 50 kpc) 
of a lO^^M© system is initially resolved by ~ 240 particles, 
and its cooling radius by ~ 6300 particles. 



2.4 Simulated observations 

In Section 3] we will examine possible discrepancies between 
the relaxed appearance of systems and the degree to which 
they are formally so. For this reason as well as to deter- 
mine the observability of features in the projected mass, 
temperature and entropy maps presented in Figs. Q to |31 
we have created mock Chandra observations of our simu- 
lations. To create these we generate 325 projected X-ray 
maps ranging from 0.5 keV to 7.0 keV in 20 eV intervals 
using the Theoretical Image Processing System (TIPSY). 
This package takes the SPH outputs from our GASOLINE 
simulations and produces smoothed projected X-ray surface 
brightness maps with the appropriate variable SPH kernel 
applied individually to the flux represented by each parti- 
cle. To these we add a three component X-ray background 
(following the prescription of the Quicksim XMM simula- 
tion package) consisting of an extragalactic power-law com- 
ponent tiChen et al. 199Z.) and two thermal components of 
10^ K a nd 10*^ K modelled with a Raymond and Smith 
plasma llRavmond fc SmithI [19771) . Generally, each image 
would then be multiplied by a energy-dependent radially 
varying effective area. However, we are interested in the de- 
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tectability of extended features at z = 0.1 and our merging 
systems rarely fit within the Chandra field of view (FOV) 
at such a redshift. For this reason, we simulate a "mo- 
saic" mode for our simulated observations using constant 
but energy dependent effective areas, area averaged over the 
Chandra FOV. An unvignetted and spectrally flat particle 
background is then added with normalisation 3.5 x 10"^ 
cts arcmin~^s~^ keV~^. The resulting images are convolved 
with an azimuthally symmetric energy dependent PSF taken 
from the on-axis PSFs specified by version 2.23 of the Chan- 
dra CALDB. The result is converted to a photon flux at a 
redshift of z = 0.1, quantised with a Poisson distribution, 
and co-added. For this paper, we use an integration time of 
50 ks. 



3 QUALITATIVE EVOLUTION 

In what follows we present a qualitative account of common 
elements in the dynamical evolution of typical cluster merg- 
ers. A series of stages which mergers generically progress 
through will be identified and the evolution of the system's 
physical and observable properties described. In the process, 
several classes of transient structures will be introduced. 
This discussion is kept brief and is intended to provide a 
context for the discussion of the observability of substruc- 
ture and relaxation which follows in Sections 2] and |Sf . For 
those seeking a more detailed description, a specific discus- 
sion of each class of transient structure can be found after- 
wards in Section |n] All quoted times are measured from to, 
the moment at which the secondary system traverses ii200 
of the primary. 

3.1 Evolutionary stages 

All of the simulations we have studied proceed through a 
similar evolutionary sequence involving five distinct stages: a 
pre-interaction phase, first core-core interaction, apocentric 
passage, secondary core accretion, and relaxation. Several 
important times during this progression which we shall refer 
to throughout our analysis are listed in Table H 

Initially, our systems are constructed to possess small 
core radii (rc ~ 50kpc, measured from /3-model fits), low 
central entropies (S — lOkeVcm^ at lOkpc) and strong cen- 
tral positive temperature gradients. The low-density outer 
atmosphere of each system becomes distorted during pre- 
interaction through tidal forces and compression, adiabati- 
cally raising the temperature of material between and pro- 
ducing a short-lived bridge in surface brightness joining 
them. As the cores continue to accelerate upon approach, 
a pair of shock fronts materialise and are driven towards 
each core, heating and compressing them briefly. The ef- 
fects of this interaction reaches its maximum strength at 
tciosest= 0.5 — 0.8 Gyr when the cores reach their closest ap- 
proach. The system's temperature and luminosity increase 
sharply for ~ 400 Myr at this point. Cooling of gas to a cold 
(T < 2.5 X 10'' K) state or to stars is subsequently quenched 

t Readers are invited to peruse relevant digital movies detailing 
the evolution of various system properties over the course of the 
simulations at http://astro.phys.uvic.ca/~babul/Merger_PaperI 



for 2 — 3 Gyr in most cases but is merely suppressed in the 
off-axis 10:1 interactions. Meanwhile, the primary system's 
projected central positive temperature gradient is strongly 
reduced and its core radius increased. A detailed analysis of 
these phenom ena is presented in the second and third papers 
of this series JPoole et alJl2006bll3. in preparation) in which 
the effects of mergers on the global properties of merging 
systems and the structure and appearance of compact cool 
cores will be studied. 

In every case (including the head-on collisions) some 
part of the secondary's cool core survives its first encounter 
with the primary core, forming a single clump and large cool 
stream of stripped material in off-axis cases and strings of 
one or more clumps in head-on cases. 

At tapo = 1.2 — 2.5 Gyr the disturbed secondary core 
reaches a maximum separation of 1.2-1.7 Mpc from the pri- 
mary core. The surviving portion then returns for another 
encounter with the primary, reaching a second pericenter 
at taccrete= 2.2 — 4.8 Gyr. In Table |5| we see that the time 
elapsed until either of these stages generally correlates with 
the mass ratio and impact parameter of the interaction, pre- 
sumably due to the declining efficiency of dynamical fric- 
tion. Second pericentric passage marks the beginning of an 
extended period of several Gyr during which material dis- 
persed from the secondary core accumulates onto the core 
of the primary system. Although a small portion of the sec- 
ondary core can survive the resulting disruption and expe- 
rience tertiary encounters with the primary, there are no 
realistically observable traces of it following taccrete in any 
of the cases we have studied^. For this reason, our systems 
are considered to be evolving as a single merger remnant 
from taccrete onwards. 

In a study of the orbital parameters of cluster substruc- 
ture in a cosmological context (including gas but not cool- 
ing), Tormon ct al. (2004) similarly find that the gaseous 
component of secondary systems becomes disrupted shortly 
after second pericentric passage. They also present simple 
analytic models for predicting the time and distance of flrst 
pericentric and apocentric passages. Our simulations com- 
pare well to their model predictions, supporting the consis- 
tency of our initial conditions with the orbital properties of 
substructure in cosmological simulations. 

Following the accretion of the secondary core, the re- 
sulting remnant proceeds through a period of relaxation. 
It appears undisturbed under reasonable observational cir- 
cumstances in simulated Chandra observations by treiax = 
4.4 — 7.8 Gyr but past this time, continues to evolve until 
the end of our simulations: recovery of the c ore's central 
temperature decrement is still proceeding (see iPoole et alJ 
■20Q,6c.., for more details) as well as accretion of both dark 
matter and gas dispersed beyond the R200 radius during the 
merger. The system sustains a virialised state within 7?5oo 
after tviriai and hydrostatic equilibrium at R500 after thydro 
(see Sectional. 

Following t closest, short-lived structures in the system's 
distribution of X-ray surface brightness, SZ signature, tem- 
perature and entropy arise. In Figs. to 1^1 we present 



' Higher impact parameters than probed by our simulations, al- 
though atypical, could result in longer lived distinct secondary 
cores. 
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projected gas surface density /X-ray brightness, tempera- 
ture/SZ and entropy maps of each of our simulations near 
several of the times discussed above. In what follows we 
shall use these figures to present the progression of these 
structures (both physical and observable) as they arise and 
dissolve. 

Readers may find it useful to refer to movies of 
our simulations* while reading the following subsec- 
tions. 



3.2 Evolution of head-on mergers 

Our head-on mergers evolve in a distinctly difi'erent way 
from our off-axis mergers, with even their qualitative evo- 
lutions exhibiting strong dependences on the mass ratio of 
the systems. We thus isolate the account of their evolution 
from that of the off-axis cases. 



3.2.1 Physical state 

In our equal mass head-on merger, two near-planar shocks 
are produced shortly before tdosest and are driven into and 
through each core (in good agreement with RSOl). They 
subsequently exit the system as two symmetric hemispher- 
ical fronts. Due to the confining geometry of this situation, 
material initially in the cores of the merging systems re- 
mains in relatively high density structures during the im- 
pact. From tciosest to approximately tapo, much of this ma- 
terial occupies an expanding disk oriented normal to the 
axis of the merging systems' motion and a filament of dense 
clumps strung between the remnants of the two cores, both 
of which are clearly visible in Fig.Q In Figs.|21andElwe see 
that this material is cold (T ~ 1.0 keV) and has low entropy 
(5* ~ lOkeVcm^). At approximately taccrete it reaccretes to 
the remnant core having radiated most of the thermal en- 
ergy generated by the impact. As a result, the system's core 
returns to a state of cooling twice as efficient at converting 
the hot ICM into stars and condensed gas {T < 2.5 x lO^K) 
as the primary prior to the merger. Although equal mass 
head-on mergers are predicted to be rare in CDM mod- 
els of st ructure formation*, this supports the findings of 
iMotl et'al . ( 2004) who find that mergers can act as a means 
of constructin g actively cooling co res. We will study this in 
more detail in lPoole et alJ (|2006cr i . 

Comparing Fig. Q to Fig. |1] we can see that reducing 
the mass of the secondary to a third of the primary's dra- 
matically changes even the qualitative behaviour of a head- 
on collision. Although two shocks are still produced, the 
one leading the secondary's motion is much more signifi- 
cant to the system's evolution. As Figs. |^ and |S| show, the 
bow-shaped geometry of this shock (again, in good agree- 
ment with RSOl) displaces the material of the primary core 
forward and laterally. The result is a conical shell of cool 
(r ~ 3keV) moderate-entropy {S ~ 150keVcm^) gas, en- 
trained by the remainder of the secondary core at its apex. 



° However, they are not without obser vational precedence; 
RXCJ0532.9-3701 jPinoguenov et al.ll2005h has boxy isophotes 
and a cloverleaf morphology in entropy, much like our relaxing 
equal mass head-on merger 



This structure surrounds a bubble of hot (T ~ 10 keV) high- 
entropy [S ~ 550keVcm^) gas. These figures also illustrate 
how the low pressure region created in the wake of the sec- 
ondary's motion is quickly and convectively filled by this 
displaced material. Material displaced most laterally reac- 
cretes first, flowing in through the back of the secondary's 
wake. Material initially displaced more forwardly is dragged 
in behind afterwards. Upon returning to the core, this ma- 
terial experiences additional shocking, resulting for example 
in the high temperature (T ~ 9.5 keV) region seen to the 
left of the core shortly after taccrete in Fig. |S| 

In the 10:1 case, the secondary core penetrates the pri- 
mary core and deposits nearly all the kinetic energy of its 
gaseous component into the core of the primary. The result 
is a large (~ 500 kpc) bubble with an expanding shell of cool 
(T ~ 3.5 keV) low-entropy (5 ~ 70keVcm^) material filled 
with hot (T ~ 7keV) moderate-entropy {S ~ ISOkeVcm^) 
gas. In Fig.lHlthe initial formation of the bubble after tdosest 
is illustrated. This bubble lasts for ~ 1.5 Gyr and is carried 
forward to a radius of ~ 750 kpc by the momentum of the 
collision. Its final extent can be seen in Figs.|S]and|U] Shortly 
before tapo it stops moving and expanding, collapses into a 
low entropy stream and accretes back to the core by taccrete. 

3.2.2 Observable state 

Following the procedure detailed in Section r2.4l we have pro- 
duced simulated Chandra observations and overlaid their 
contours onto the projected gas surface density maps pre- 
sented in Figs. Qll] and This allows us to determine which 
distributed features have sufficiently high surface bright- 
nesses to be detected under reasonable observational circum- 
stances. We have also generated SZ maps and overlaid their 
contours on the temperature maps presented in Figs. |21 El 
andlHI Since these maps do not include the effects of instru- 
mentation, we cannot comment in detail on the detectability 
of low surface brightness features through SZ observations. 
We place our faintest co ntours at a level of lo g y — —5.5 
{5T = -17^^" at 30GHz. iMcCarthv et alJl2003t) to approx- 
imate the depth to which the next generation of SZ imag- 
ing instruments w ill reach {e.g. —lO/iK for SZA at 30GHz, 
iKnox et al lEooi). The resulting SZ limit is a good match 
to the system extent realistically observable by Chandra. 

In the 1:1 case, we can see in Figs. and |5| that both 
the expanding disk and filament of disturbed core mate- 
rial created in the impact and discussed above are rea- 
sonably observable through X-ray observations. Comparing 
these figures we can see that the remnant's X-ray isophotes 
are extremely boxy while the SZ contours are significantly 
more regular and much less elliptical (this is a trend that 
we observe in all of our simulations). Complicated temper- 
ature and entropy variations are seen throughout the du- 
ration of the merger, with small ~ 100 kpc scale variations 
(AT/T ^ 15%) remaining at treiax- 

In the 3:1 head-on case. Figs. ^ and illustrate that 
the extended features generated in this merger following 
tcioseat are similarly observable as in the 1:1 case. However, 
we can see that at taccrete there are large extended low sur- 
face brightness features that escape detection by Chandra. 
Complicated temperature and entropy variations, distinctly 
different from those observed in the 1:1 case, emerge after 
taccrete and pcrslst Until the end of the simulation. The tem- 
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perature fluctuations present at t^eiax are more significant 
[AT/T ^ 20%) and spatially irregular than in the 1:1 case. 

Lastly, Figs. |7| and |S| illustrate the observability of ex- 
tended structures formed in our 10:1 head-on merger. We 
can see from these that the large expanding bubble formed 
from the collision is not detectable when it reaches its largest 
extent at tapo- Furthermore, the system looks quite regular 
in both X-ray and SZ at t accrete, although small but signif- 
icant irregularities in the central X-ray isophotes persist for 
~ 2 Gyr afterwards. Temperature and entropy fluctuations 
are present throughout the merger with the bubble formed 
from the core material dominating the temperature and en- 
tropy distributions after tdosest- Significant fluctuations in 
temperature {AT/T ~ 20%) are present at treiax, with little 
change occurring between taccretc and trciax- 

3.3 Evolution of off-axis mergers 

In all of the off-axis simulations we have studied, several 
distinct classes of transient structure are formed when the 
secondary (and sometimes primary) system's core becomes 
disturbed. Interestingly, the qualitative evolution of these 
structures is remarkably independent of impact parameter 
and mass ratio over the ranges we have studied. Further- 
more, their properties are consistent with those of "cold 
fronts" observed in man y Chandra and XMM observations 
jMarkevitch et all200(il) . Cold fronts are contact discontinu- 
ities between regions of bright, dense, cold gas and faint, hot, 
rarefied gas. These regions are observed to be in near pres- 
sure equilibrium with the surrounding ICM and the jumps 
in temperature across their interface imply mildly sonic or 
sub-sonic velocities. 

Several authors have studied simulations of merging 
clusters illustrating that the ram pressure disruption of cold 
low entropy cores during mer gers can account f or these 
featu res (Bialck ct al. 2002; Hei nz et al] l2003: Mat his et alJ 
[200^. Once removed from the dark matter potential confin- 
ing it, stripped gas lags behind its dark matter and galac- 
tic components, adiabatically expanding and cooling until it 
establishes pressure equilibrium within the resulting merger 
remnant. 

Our simulations confirm these results. In what follows 
we will give a qualitative account of our off-axis mergers' 
evolution through a series of transient structures and dis- 
cuss their observable properties. We shall find that there 
are many points when several classes of structure resembling 
cold fronts can be created. Additional details can be found 
in Section |S] where each of these classes will be studied in 
detail and new mechanisms involved in the production and 
evolution of cold fronts will be presented. 

3.3.1 Physical state 

In all of our 1:1 and 3:1 off-axis mergers, ram pressure strip- 
ping causes both the primary and secondary cores to develop 
"comet-like" head/tail morphologies in gas surface density. 
X-ray surface brightness, temperature and entropy during 
their first pericentric passage. In the 10:1 cases, this is true 
only for the secondary core. In all cases, this morphology 
starts ~ 100 Myr before tdoseat and ends shortly before tapo 
with a typical duration of ~ 0.5 Gyr. 



Upon reaching apocentric passage, several extended 
structures in the system's distribution of temperature, en- 
tropy and gas density have evolved. In Figs. Q and 2] we 
see that a bridge connecting the surviving portions of the 
primary and secondary cores has formed in the 1:1 and 
3:1 cases (see Section \<6.1l . Several long and roughly radial 
outwardly moving plumes of cool (T ~ 3keV) low-entropy 
(S ~ GOkeVcm^) material have also formed (see Section 
10.21 . These features are not obvious in the system's gas den- 
sity or surface brightness maps, but are evident in the tem- 
perature and entropy maps presented in Figs. and |S| 
All of these structures are essentially absent in our off-axis 
10:1 mergers. 

Following tapo, material displaced as plumes from the 
primary and secondary cores stops expanding away from the 
system and begins returning towards the remnant core. At 
taccrete, the Secondary core returns for its second (and last) 
pericentric passage of the system. After this, plume mate- 
rial begins accreting to the remnant core as coUimated high 
velocity (~ 1000 kms~^) streams of material (see Section 
10.41 . At this point, evidence of a secondary core has been 
destroyed and these features are not associated with an ap- 
parent merging system. Stream material is generally shock 
heated as it accretes to the remnant core. These streams 
occur even in our 10:1 cases and are visible in the tempera- 
ture and entropy maps presented in Figs. and |5| 
as narrow and radially oriented cold low entropy structures 
following taccrete. Wheu they accrete directly to the rem- 
nant core, these streams lead to multiple discontinuities, or 
"edges" in the system's central surface brightness distribu- 
tion (see Section fO.Sl . 

3.3.2 Observable state 

Elongated and curved disturbed cores have been observed 
in several systems (e.g. A85, A1758, and A2256) and are 
frequently attributed to the disruption of accreting systems 
due to ram pressure stripping. We observe this for only a 
short duration (< 0.5 Gyr) following tdosest in all of our 
off-axis mergers. Prior to tdosest, this morphology does not 
manifest. At tapo, longer lived elongated tails are created 
but the ejection mechanism is not ram pressure stripping 
(see Section [0.21 . 

In Fig. we see that the bridge formed after tdosest 
in our off-axis 1:1 mergers remains visible in our simulated 
Chandra observations until destroyed at taccrete. The mate- 
rial forming this structure is not stripped from the central 
core but from the region just outside the initial cooling ra- 
dius. Hence, it is of intermediate entropy (300 — 500 keVcm^) 
and temperature (~ 6.5 keV). No strong pressure gradients 
arise across it, as shown in the SZ maps in Fig.|21 

The bridges formed in our off-axis 3:1 mergers remain 
detectable for only 1.0 Gyr following tdosest, or roughly until 
tapo- As in the 1:1 cases, the entropy of the bridge material is 
of intermediate entropy (200—400 keVcm^) and temperature 
(~ 4.5 keV) and exhibits no strong pressure gradients. 

In Figs. Q and |1] we also see that the plumes generated 
from the disturbance of the primary and secondary cores are 
observable in the 1:1 cases but only marginally so in the 3:1 
cases. Figure Q shows that in the 10:1 cases, although the 
system appears disturbed after tdosest, no extended features 
exist with sufficient surface brightness to be detected. 
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Figure 10. Displacement of the centre of mass of gas within the central 150 kpc from its dark matter centre of mass (smoother, 
lower-amplitude red curve) and the mass-averaged velocity of the same material in the direction of that displacement (the more erratic 
blue curve) for the late stages of each of our simulations. Vertical dashed lines correspond to taccrete and t^^i^^. The increased activity 
between these two timescales is primarily the result of the gas core being disturbed by the secondary's gas core during the second 
and typical the final pericentric passage, and by the ensuing cold stream. Values across the top indicate the transverse velocity at to 
depicted by each column while text on the right indicates the mass ratio depicted by each row^ . 



At taccrete, Very Significant temperature fluctuations 
(AT/r ~ 20%) are present in all cases and persist until 
well after treiax- Because the accreting streams retain angu- 
lar momentum from the merger and do not readily mix with 
the remnant, they produce spiral structures in the remnant's 
temperature and entropy distribution which are visible in 
our temperature and entropy maps at treiax- 

Because the SZ signature is proportional to gas den- 
sity (pg) rather than the Pg dependence of the X-ray surface 
brightness, extended low density features tend to be washed- 
out in the SZ maps. Despite this, we can see that our SZ 
maps indicate that the bridges formed in the 1:1 and 3:1 
mergers remain above our nominal SZ threshold. At t,.eiax , 
when accreting streams leave the X-ray surface brightness 
elongated in these cases, the SZ maps look significantly more 
regular with only the 1:1 Vt/Vc = 0.15 case looking compa- 
rably disturbed to its X-ray surface brightness. 



3.4 Core oscillations 

In several of the cases we have studied, we find that the 
central (r < Q.S-Rsoo) ICM of our merger remnants oscil- 
late with respect to their dark matter distributions. In the 
discussion of the evolution of disturbed morphology and re- 
laxation which follows, there will be several instances when 
these oscillations will be relevant. Such oscillations have 
been discussed by several authors in attempts to account for 



"edges" observed in the central surface brightness of some 
clusters (see Section fti. 51 for more details and references). 

In Fig. llOl we present the displacement of the centre of 
mass of the densest central 150 kpc of our remnant gaseous 
cores relative to the centre of mass of the densest central 
150 kpc of dark matter. To differentiate oscillating modes 
from the infiuence of substructure passing through the core 
in this calculation, we also present the mass averaged veloc- 
ity of the gas in the direction of this displacement's vector. 
Displacement due to coherent oscillations will be accompa- 
nied by positive velocities until peak displacement when the 
velocity turns around and becomes negative, returning to 
zero when the displacement does. The displacement of the 
core due to the arrival of substructure is revealed by in- 
creases in the displacement at times of small or negative ra- 
dial velocity. There are several instances of significant core 
oscillations following taccrete stimulated by the second peri- 
centric passage of the secondary core. They generally last 

until treiax- 

We have performed this calculation for material within 
50 kpc as well. On these scales we observe no oscillations 
larger than the softening length (10 kpc) at times when os- 
cillations are present on 150 kpc scales. Hence, gas in the 
central regions of our remnant cores (where cooling is most 
rapidly taking place) remains tightly coupled to the dark 
matter potential with the surrounding material constituting 
the oscillation seen on 150 kpc scales. It has been hypoth- 
esised that energy from bulk oscillations in the potential 
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of the core due to displacements of the gas from the dark 
matte r could act as a source of heating for co mpact cool 
cores fMarke vitch. Vikhlinin. fc M azzottal l200 J) . Our anal- 
ysis suggests that although oscillations may occur on 150 kpc 
scales, the actively cooling gas on 50 kpc scales can remain 
strongly coupled to the dark matter. In such cases, oscilla- 
tions in the potential may not be affecting cooling. Further 
study is required to determine if this is the case for oscilla- 
tions stimulated through close encounters by secondary sys- 
tems with mass ratios of only a few p ercent, such as those 
studied bv lXittlev fc HenriksenI J2005l) . 

All of these calculations determine the motion of the 
gaseous core relative to the dark matter core. We have 
checked to see how the dark matter core is itself moving 
within its halo. For our 1:1 cases, we find that on both 150 
and 50 kpc scales the dark matter cores remain stationary 
relative to the centre of mass of the system following taccrete- 
However, the dark matter cores of our 3:1 and 10:1 mergers 
experience significant motions about the system's centre of 
mass until well after treiax- The oscillations we observe be- 
tween the gas and dark matter on 150 kpc scales are thus a 
result of the dark matter core orbiting the centre of mass 
of the system. As the gravitationally dominant dark mat- 
ter core moves, the gaseous component is dragged behind 
it but with a lag due to the extra pressure it experiences. 
The result is an oscillation between the two components as 
the dark matter core moves in its orbit. On 50 kpc scales, 
baryonic material dominates the potential of most of the gas 
and the dark matter can follow it since it is not subject to 
pressure. As a result, oscillations between the gas and dark 
matter do not occur on these scales. 

Following taccrete, the gascous and dark matter cores 
of our 1:1 merger remnants rapidly become concentric by 
treiax- A fcw displacements of the core position with respect 
to the dark matter are apparent but these are generally ac- 
companied by negative radial velocities indicating that they 
are a product of the accretion of substructure onto the core. 

The cores of our 3:1 merger remnants all experience 
bulk oscillations between their gaseous and dark matter 
components following taccrete- In the head-on and vt/Vc ~ 
0.15 cases, the amplitude of these motions are 100 kpc and 
150 kpc respectively with velocities reaching 400 kms~^. 
These motions last for ~ 1 Gyr and end immediately before 
treiax- lu the Vt/Vc ~ 0.4 case, a series of 4 core oscilla- 
tions occur beginning at taccrete with amplitudes as large as 
50 kpc and velocities of 175 — 250 km s~^ . They are of much 
shorter duration (~ 0.5 Gyr) and are finished by treiax- 

Some part of the dark matter cores of the secondary 
systems in our 10:1 mergers remain intact in near circular 
orbits within 50 kpc of the remnant's centre until the end 
of our simulations. This is not surprising since the sinking 
timescale due to dynamical friction is the longest. As a re- 
sult, we measure small amplitude (< 20 kpc) oscillations in 
our 10:1 merger remnants. These motions are most coher- 
ent in the head-on case and involve bulk velocities of < 100 
kms~^. 



4 EVOLUTION OF DISTURBED 
MORPHOLOGY 

With the general evolution of merging systems accounted 
for, we now discuss the timescales by which evidence of 
the merger event disappears. Determination of the degree 
to which a cluster is disturbed has generally taken three 
approaches in observational studies: visual inspection, mea- 
surement of displacements of the X-ray peak from the X-ray 
or mass centroid of the system, and power-ratio analyses 
of the system's X-ray surface brightness. In this section we 
shall apply each of these approaches in turn to establish 
when our simulated merger remnants would return to ap- 
parently undisturbed states. In the section which follows, 
we shall compare these timescales to more formal and phys- 
ical measures of the system's dynamical state. 

4.1 Simulated observations 

In evaluating the effectiveness of quantitative measures of 
our remnants' dynamical state, we shall use as a fiducial the 
time when our systems appear relaxed under reasonable ob- 
servational circumstances. We determine this time by visual 
judgement of our simulated Chandra observations and of the 
uniformity of their isophotes; a process we facilitate by gen- 
erating unsharp masked images with smoothing kernel sizes 
of 1.5" and 15" (Chatzikos ct al. 2006). This is all done in 
projections normal to the plane of the orbit. Substructure 
is likely most identifiable in this projection and thus, this 
measure represents a conservative upper limit on the time 
when the system would look disturbed. We have studied how 
our results change if evaluation is performed in a;-projections 
and find that the systems tend to look undisturbed as much 
as 2 Gyr earlier. 

Generally speaking, evidence of substructure disappears 
first from the simulated images, then from the unsharp 
masked images and finally in the structure of the X-ray 
isophotes. The delay between the dissolving of substruc- 
ture in the images and the regularisation of their isophotes 
can take several Gyr and varies significantly from case to 
case. Since the isophotes are the most sensitive measure, 
the system will be said to appear relaxed at treiax when 
all isophotes brighter than 6x the background level (ie. all 
but the faintest in Fig.Q^JS'iidlTJ appear regular by visual 
inspection. These times are tabulated in Table |5| 

In one case (1:1 head-on), minor complications arose 
during this procedure. Although the simulated observations 
appear regular and relaxed hy t ^ 3.8 Gyr, their isophotes 
remain significantly boxy until t ~ 4.4 Gyr. Some evidence of 
this boxy morphology lingers until the end of the simulation 
making a clear determination of treiax difficult. 

4.2 Centroid shifts 

The apparent displacement of a system's core from the cen- 
tre of the system has been used as a measure of dynami- 
cal disequilibrium by several authors and has been imple- 
mented in several different ways. For instance, Mohr et aP 
il993h quantify X-ray centroid shifts (wx) using the offset 
of a system's X-ray surface brightness peak from its sur- 
face brightness centroid (integrated within a projected ra- 
dius rp). They argue that systems for which Wx is dependent 
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Figure 11. Evolution of the centroid variance computed from the X-ray surface brightness maps in the x (line of sight) and z 
(orthogonal) projections (thin re d lines and thick blue lines respectively). Horizontal dashed lines indicate the range of centroid shifts 
observed bv lO'Hara et all i200(i) while vertical dotted lines indicate (from left to right) to, t closest i ^accrete and t^^i^^. Centroids 
are computed within 1.2 Mpc (roughly the initial -R500 of the primary). The centroid statistic traces the qualitative evolution of the 
merging system very well. The statistic takes on its maximum values when the system is in its most disturbed state and typically, it 
falls off b etween tgccrete and trsJg.T as the system recovers. Over the course of the merger, the statistic spans the entire range of values 
found bv lO'Hara et all l200ffl in an analysis of 45 clusters^. 



on Tp are dynamically young. For a sample of 46 Einstein ob- 
served clusters, they obtain an average value (wx) ~ 62kpc 
with 71% of systems having > s.t a, 3a or greater sig- 
nificance. _ITiese_j;esults_gener with more recent 
results of iKolokotronis et alJ l|200!j) who find (wx) = 83 kpc 
for a sample of 22 ROSAT observed systems. 

We have chosen to quantify the ce ntroid variance of ou r 
systems using the method employed bv lO'Hara et alJ i2006l) . 
This approach determines the offset between the X-ray peak 
and centroid as a function of < -R500 and assigns to Wx 
the RMS of the resulting profile. We determine this profile 
by computing the centroid within a circular aperture of size 
rp ~ -R500 and then reducing its size by 5% until it reaches 
rp ^ O.OS-Rsoo- We have experimented with various details 
of this procedure and find that it is essential to exclude 
the bright central core if the statistic is to detect subtle 
distortions. Hence, we excise the central rp — 30 kpc when 
computing centroids (but retain it when we compute the 
position of the peak). 

In Fig. Illl we plot the time dependence of Wx measured 
in this way for each of our simulations in x and z-projections. 
The evolution of this quantity in y-projections is very similar 
to its evolution in the 2-projection and is omitted to clarify 
the plot. We also indicate th e range of typical values ob- 
served bv lO'Hara et all (l2006l) in their study of 45 ROSAT 
observed clusters {wx = O.OO8-R500 to O.OT-Rsoo, with 5 ex- 



treme points excluded) as well as the times of several im- 
portant stages in each system's evolution. 

As a product of our idealised initial conditions, all of our 
systems start with values of Wx less than the lower bound of 
the observations. In off-axis cases, it immediately begins to 
increase in a;-projections as the cores slowly move apart in 
projection. This effect grows more significant with increas- 
ing Vt/Vc as we would expect. In a-projections, no increase 
occurs until ~ to when the body of the secondary enters the 
R500 ~ 1.2 Mpc aperture used for our calculations. In the 
off-axis cases, the statistic for the x and z-projections are 
virtually indistinguishable following tdosest- In the head-on 
cases, the signal is much lower when the merger proceeds 
along the line of sight. Typically, the statistic drops sharply 
when the primary and the secondary are close together, as in 
the line-of-sight/head-on case and more generally at tdosest 
and t accrete- This is due to the fact that the centroid shift 
in all the apertures is nearly the same. In detecting distur- 
bances, the centroid variance method depends on the dis- 
placement varying as the aperture is reduced and succes- 
sively more distant surface brightness fluctuations are ex- 
cluded from consideration. 

Our 1:1 mergers evolve to values of Wx much larger than 
those observed for a duration of ~ 4Gyr. During this time, 
the system is a distinct binary a nd is not repre s entati ve of 
any of the systems present in the lO'Hara et al"] feOOd) cat- 
alogue. Following taccrete, Wx decliues rapidly and reaches 
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Figure 12. Time evolution of the displacement between the X-ray peak and the centroid of the projected mass maps in the x (line 
of sight) and z (orthogonal) projections (t hin red lines and t hick blue lines respectively). Horizontal dashed lines indicate the mean 
and maximum centroid shifts measured bv lSmith et al.l J2005ll (normalised by -R500 of our system) while vertical dotted lines indicate 
(from left to right) to, t^-iasest i taccrete and t^elax- In the case of massive (1:1) mergers, the statistic behaves like Wx and tracks the 
morphological evolution of the system over the course of the merger. For moderate and minor mergers, the statistic remains high well 
past t^^iax although the system is close to being virialised for t > t^siax- This signal is entirely due to low-mass remnant dark matter 
substructure that keeps jiggling the mass centroidt. 



a roughly constant level near the minimum observed value 
~ 1 Gyr before treiax- In our vt/Vc = 0.4 case, there is a sub- 
sequent rise at late times due to a late accretion of material. 
It is likely that the strength of this rise would not be as 
significant in a system hosting realistic substructure where 
the effect of late accreting material from multiple mergers 
would tend to cancel their individual effects. For this rea- 
son, we will consider the system to have relaxed according 
to this measure ~ 1 Gyr prior to treiax in this case but ac- 
knowledge that more study is required to understand clearly 
the role long lasting substructure has on Wx in realistic en- 
vironments. 

In the 3:1 cases we see a similar rise after to to val- 
ues greater than the maximum observed but lasting for a 
much shorter duration ('^ 2 Gyr). Following taccrete, Wx os- 
cillates erratic ally a mongst the highest values observed by 
lO'Hara et alJ i2006l) . It slowly declines to the minimum ob- 
served value by the end of the simulations in the head-on 
and Vt/Vc = 0.4 cases but remains high in the vt/Vc = 0.15 
case. The cause of this slowed or absent decline in Wx follow- 
ing taccrete iS a sustained shift of the centroid in the central 
0.5i?5oo- This is due to oscillations in the remnant's central 
dark matter distribution. 

Lastly, our 10:1 cases are disturbed to high levels of Wx 
for ~ 4 Gyr following to and for ~ 1 Gyr following taccrete 
in the Vt/Vc = 0.4 case (due to the reduced disruption of 



the core). They oscillate just above the minimum observed 
value otherwise. 

Thus, our systems appear exceptionally disturbed to 
our centroid variance measure for ~ 4 Gyr between tdosest 
and taccrete wheu they generally appear as two distinct clus- 
ters. Immediately following taccrete for our major/moderate 
(1:1 and 3:1) mergers and ~ 2 Gyr after tdosest for our 10:1 
mergers (when the systems appear as single remnants), our 
simulations all exhibit values of Wx similar to the maximum 
observed by O'Hara et al. (2006) (with the natural excep- 
tion of the head-on cases in a;-projections). In 3:1 cases, Wx 
remains high until well past treiax - Our equal mass and 10:1 
mergers quickly rel ax to levels compara ble to the minimum 
values observed bv lO'Hara et alJ ll2006l) . Referring to Table 
|5|we can see that for a cluster observed at z = 0.1, most 
10:1 mergers occurring before z ~ 0.5 will have evolved past 
taccrete and wiU thus yield Wx ^ O.Ol-Rsoo- This fact com- 
bined with the expectation that most clusters should have 
experienced such an event naturally explains the observed 
minimum of this statistic. 

We have also computed this statistic for SZ maps {wsz) 
of our systems. We find that the evolution of wsz is nearly 
identical to that of Wx, with both producing very similar 
maximum shifts and periods during which Wx is significantly 
higher than O.l-Rsoo- 

Another ap proach to measuring centroid shifts has been 
implemented bv lSmith et alj (|20o5 who measure the shift 
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between the peak in the X-ray surface brightness and the 
centre of projected mass of the system (Ar), as determined 
from weak lensing maps. For a sample of 10 clusters, they 
find (Ar) = 55kpc (~ 0.045-R5oo for our systems ) with 70% 
of their systems having Ar > at a 3cr or greater signifi- 
cance. In Fig. 1121 we present the evolution of this statistic 
for our simulations. Again, we plot only x and ^-projections 
due to the similarity of results in y and z-projections. 

Examining Fig. 1121 we see that for the 1:1 mergers, 
this statistic drops precipitously after taccrete (similarly to 
Wx) and sustains similar values after treiax- For the 3:1 
mergers however, there is a long lived late time shift in z- 
projections which does not reach maximum amplitude until 
treiax and does not reduce to Ar ~ 0.01i?5oo until shortly 
before the end of our simulations. A similar but even longer 
lived shift persists in our 10:1 mergers which fail to sus- 
tain Ar <^ O.OI-R500 for more than 0.5 Gyr in any projection 
for any case. These sustained shifts are due to central dark 
matter oscillations (see Section [3.41 . 

Given the relatively high frequency of 10:1 mergers and 
the long duration of their impact on Ar, we expect that few 
systems should exhibit Ar < 0.045ii5oo (the mean value for 
our 10:1 vt/Vc ~ 0. 15 merger af t er tre iax) and therefore, 
it is surprising that ISmith et al.l ll2005l) find any systems 
consistent with Ar = 0. However, if projection effects are 
considered, it may be that some or all of the 3 systems to 
which they assign relaxed morphologies using this statistic 
are being seen with their X-ray peak merely in projection 
against their centre of mass. Alternatively, long lived sub- 
structure from multiple past minor mergers may compete 
against each other in their infiuence on the centre of mass, 
reducing its overall offset from the X-ray peak. Further study 
of this statistic and its efficacy as an indicator of disturbed 
systems in a full cosmological context is clearly necessary. 

4.3 Power ratio morphology 

The use of surface brightness moment decomposition for 
quantifying the degree to which sys t ems appear relaxed 
was first introduced bv lBuote fc Tsaa jlOOj) an d has since 
been utilised in sev eral observational llBuote fc Tsai 1 996; 
Ijoltcma ot al. 2005') and theoretical studies jThomas et all 
Il99a: .O'Hara et al..,206i) . This method uses power ratios 
(Pm/Po) computed from moments of the system's surface 
brightness distribution; E(a:), where x' = ( R, 4>) is in po- 
lar co ordinates. Following the procedure of Ijeltema et alJ 
i2005l J05 hereafter), we have applied this procedure to our 
simulations computing power ratios from 

Po{R) = [a„\n{R)f (2) 

P^(R) = l^J^i^l^^^'r^) (3) 
where and bm are moments of the distribution given by 

am{R)= [ T.{x')(R')"' cos m(l>' (fx' (4) 

J R'^R 

bm[R)= ( T.{x'){R')"' smm(j}' (fx' (5) 
Jr's^r 



These statistics are usually normalised by Po («-e. expressed 
as Pm/Po), a measure of the total luminosity of the system. 

We have studied the effects of removing the central cores 
of our simulated clusters when computing these statistics. 
We find this to be particularly necessary when the system 
is only subtly disturbed. At such times, the contribution of 
any persisting structure in T.{x) to Eqns. 4 and 5 becomes 
overwhelmed by the bright and centrally concentrated cores 
of the merging systems. When we exclude the central cores 
from our calculations of Pm we find a significant increase 
in amplitude and scatter. For the discussion which follows, 
we have computed power ratios with the central 30kpc ~ 
0.025i?5oo excised. We have confirmed that there is little 
effect on our conclusions if this radius is reduced to 15 kpc 
or increased to 50 kpc. 

In Fig. 1131 we plot the time evolution of P2/P0, P3/P0 
and P4/ P0 in a 0.5 Mpc z-pro jected aperture for each of our 
mergers. iBuote fc Tsail il995ll find that model systems with 
power ratios identically equal to zero can be assigned signif- 
icant values (a few times 10~^ for P3/P0 for instance) when 
realistic noise is added to simulated X-ray images. Based on 
their analysis of an observed sample of clusters, J05 found 
that median noise contributions resulted in values for P2/P0, 
P3/P0 and P4/P0 of 7.5 x 10"®, 1.9 x 10"*, and 9.4 x 10"'' 
respectively. We indicate these levels in the plot. The plot il- 
lustrates that in the absence of noise, the values of the power 
ratios for the mergers are typically well below the median 
noise levels, except when the systems are in their most dis- 
turbed state. In practise, the values assigned to a relatively 
relaxed-looking cluster will be highly dependent on the qual- 
ity of the individual observation. In fact, the lowest values of 
Pn / Po quoted by J05 are a factor of ~ 10 below the median 
noise levels. In what follows we shall loosely refer to the me- 
dian noise levels reported by J05 as the "observable levels" 
of each statistic. We will treat these levels as a reasonable 
fiducial when determining the point at which our systems 
would be judged as undisturbed by this approach. 

The evolution of P3/P0 and P4/P0 is very similar; the 
latter behaves more or less like a scaled version of the former. 
In the approach to the first pericentric passage at tdosest, 
both power ratios experience a brief sharp spike which ex- 
ceeds observable levels in all cases, though just barely in 
the 3:1 and 10:1 mergers. This is followed by a climb to a 
plateau or a broad peak, during which time the signal re- 
mains at or below the observable levels except in the 1:1 case, 
and a fall-off at ~ taccrete ± 1 Gyr. The amplitudes of the 
initial spike and of the plateau/broad peak are correlated 
with the secondary's mass, with the most massive merger 
generating the greatest signal. There does not appear to be 
any clearly discernible correlations with impact parameter. 
Once the signal falls below observable levels after the second 
drop, it remains so for the duration of the simulation under 
all conditions. 

The evolution of P2/P0 is different from that of P3/P0 
and Pi/Po- Like the latter two, the signal brieffy spikes 
above the noise level at tdosest but rapidly drops below. 
It exceeds the noise level again at taccrete for < 2 Gyr and 
then generally only for the mergers with mass ratios greater 
than 3:1. (At 3:1, the second spike just exceeds the observ- 
able level). Afterwards, it stabilises at values far below ob- 
servable levels, changing very little after treiax- For the ma- 
jor/moderate (1:1 and 3:1) mergers, P2/P0 settles at ~ 10~^ 
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Figure 13. The time evolution of the X-ray surface brightness power ratios computed within a 0.5 Mpc aperture in the z-projection, 
the projection in which the merger-induc ed disturbance s are t he most apparent. Dotted horizontal lines depict the median noise levels 
for the observational sample analysed bv lJeltema et af ]l|2003) for (in descending order) P2/P0, P3/-P0 and Pi/Po. Dash-dotted black, 
thick red and blue lines correspond to P2/PQ, P3/P0 ^'Hd Pi/Po respectively. Vertical dotted lines represent (from left to right) to, 
tcloseat! taccrete and t^f,iax- The power ratios can easily pick out the very short-duration highly disturbed states but for the most part 
the amplitudes tend to be comparable to or less than the median noise levels for the J05 sample of clusters^ . See also Fig. 1141 



while our 10:1 mergers, the power ratio hovers at values of 
about ~ 10~^. 

In their study of the redshift evolution of power-ratios, 
JOS find that between z = 0.5 — 0.9 (their high-2 sample) and 
z < 0.5 (their low-z sample), the P3/P0 and Pi/Po ratios 
show statistically significant declines while the P2/P0 ratio 
shows no significant change. Furthermore, observed values 
of P2/P0 are confined to a range more narrow than P3/P0 
or Pi/Po. Our simulations reproduce these trends but with 
a much lower normalisation. Further study is required to un- 
derstand if this discrepancy is a product of noise effects or 
our idealised initial conditions and lack of cosmological sub- 
structure. Since P2/P0 is the power ratio most sensitive to 
the large scale ellipticity of the system, it would be interest- 
ing to determine if more realistic triaxial initial conditions 
would produce a similar narrow range of near constant final 
P2/P0 but with a higher normalisation comparable to that 
observed. As it stands, our results suggest that the power 
ratios are most sensitive to gross disturbances such as those 
that arise early in the merging process. 

Correlations between different power ratios have gener- 
ally been used to identify disturbed systems and to study 
the evolution of cluster substructure. Here, we focus our at- 
tention on the P4/P0-P3/P0 plane. In Fig. 1141 we present 
a comparison in this plane between the observations of J05 
and our simulations computed within apertures of 0.5 Mpc. 
Results computed in x, y and z-projections are presented 
with red triangles, green squares and black circles respec- 



tively. In Fig. 1141 we can see that the observed trends in nor- 
malisation, maximum amplitude and scatter are well repro- 
duced, though our measurements cluster at values of P3/P0 
and Pi/Po well below the observations. These two power 
ratios return values above the observable levels only during 
the highly disturbed phases of the merger, and during the 
initial stages of the off-axis mergers. 

We have also examined the distribution of points in the 
P4/P0-P3/P0 plane, where the power ratios are computed 
within a larger aperture of radius -R500 1.2 Mpc. With the 
exception of some reduction in the scatter, the results are 
very similar to those presented in Fig. 1141 

Finally, we also applied the power ratio analysis to total 
projected mass and SZ maps to examine the efficacy of this 
approach to measuring substructure through weak lensing 
and microwave observations. The results were indistinguish- 
able from the X-ray results, with amplitudes of the P3/P0 
and Pi/Po for the projected mass maps being lower. 



4.4 Comparison of disturbed morphology 
measures 

We have examined several methods of determining the dy- 
namical state of galaxy clusters. In Figs. [TTI 1121 and 1131 we 
present the time evolution of Wx, Ar and Pm/Po- In each 
case we indicate the moment when the isophotes of our sim- 
ulated 50ks 2 = 0.1 Chandra observations appear undis- 
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Figure 14. Plot of P4/P0 vs. P3 / Pq for our simulations measured at different time snapshots over the course of the merger, within an 
aperture of 0.5 Mpc. The data of Ijeltema et al. (2005) is plotted in blue with our simulations plotted as red triangles, green squares 
and black circles for x, y and 2 projections respectively. This combination of power ratios yields a signal above the median noise level 
of JOS sample (shown as horizontal and vertical lines) only when system is caught in the short-lived highly disturbed state. It is most 
sensitive to disturbances induced by highly off-axis mergers and least sensitive to line-of-sight mergers^ . 



turbed as a fiducial reference for the sensitivity of each ap- 
proach. 

We find that while the centroid shift and power ratio 
approaches are both capable of identifying highly disturbed 
systems or systems with significant, well-defined substruc- 
tures, the Wx is more sensitive to subtler disturbances and 
its time evolution tracks the visual appearance as well as the 
dynamical state of the merging system quite well. The am- 
plitude of the power ratios is relatively low when the system 
is not highly perturbed (for example, the analysis barely 
registers 10:1 mergers) and appears to be highly suscepti- 
ble to noise in the X-ray images. More study is required to 
understand in exactly what way realistic noise affects the 
sensitivity of the Pm/Po statistics. 

Our experime nts with the Ar statistic introduced by 
ISmith et al.l lj200^ suggest that it is very sensitive to long 
lived substructures and oscillating modes in the dark mat- 
ter haloes of merger remnants. The statistic yields a signal 
even when dynamically the system is within 2% of being 
relaxed as defined by our virial indicator (see Section |K| and 
Fig. I15|l . However, it is likely that our idealised mergers are 
exaggerating this sensitivity. More study with cosmological 
simulations which include a proper account of halo substruc- 
ture should be conducted to determine exactly how effective 
Ar is relative to centroid shift and power ratio measures. 

In summary, for clusters ai z — 0.1, the eye is remark- 
ably adept at identifying disturbed morphologies in deep 
Chandra images but quantitatively, we find that the cen- 



troid variance, w^, is the most effective measure of a sys- 
tem's dynamical state. 



5 RELAXATION 

Is a system which appears undisturbed actually in a relaxed 
state? This is commonly assumed to be the case in studies 
which have selected compact cool core systems with reg- 
ular and symmetric isophotes to represent equilibrium sys- 
tems. However, all clusters are a product of mergers and even 
though their dynamical timescales are short (~ 1.5 Gyr at 
R200), the timescales by which bound substructures become 
disrupted and dissolve into the accreting halo can be much 
longer, as illustrated in the preceding sections. As a result, 
a merger remnant can take significant lengths of time to 
reach a proper equilibrium state, perhaps even longer than 
the typical interval separating merger events. 

In practical terms, the point at which a system appears 
undisturbed depends on the state of the system's baryonic 
component whose signatures of disturbance can be diluted 
by observational limitations. Furthermore, the gaseous com- 
ponent is subject to disruptive gas dynamical forces which 
do not act on the dynamically dominant dark matter com- 
ponent, perhaps erasing evidence of substructure faster than 
it loses dynamical relevance. Hence, it is natural to suspect 
a discrepancy between when a system appears relaxed and 
when it is formally so. 

In this section we will quantify the relaxation of our 
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merger remnants by studying their recovery towards a viri- 
alised and hydrostatic state of equihbrium. We shaU com- 
pare the timescales by which they do this to the point at 
which the system would look relaxed under typical observa- 
tional circumstances. We shall find that our systems typi- 
cally achieve a virialised state approximately when they ap- 
pear relaxed but continue to exhibit deviations from hydro- 
static equilibrium at the level of 10 — 20% until the end of 
our simulations. 



5.1 Virialisation 

To study the virialisation o f our systems we employ 
the scalar virial theorem (see ICoUinsI Il978l : ISoitzerl Il978t 
iBinnev fc Tremaind HqStI for good reviews) which states 
that 



2 dt2 



= 2T + W + S 



(6) 



where T = K -\-U is the sum of the system's kinetic [K) and 
thermal {U) energies, W = Xiili Fi ■ fl is, the 'virial', and 
/ is the system's moment of inertia, all computed within a 
specified volume. The term denoted hy S = Sg + Sd is & sur- 
face pressure term integrated over the bounding surface of 
the volume, with contributions Sg and Sd due to the gaseous 
and coUisionless (dark matter and stellar) components given 
by 



Sg ^ - i PTdA 



Sd = 



(7) 
(8) 



where P = nkT is the pressure of the gas with a number 
density n and temperature T and <jd and pd denote the 
velocity dispersion normal to the surface and mass density 
of the coUisionless component. 

Averaged over an interval in which the system is peri- 
odic (or instantaneously for a steady state system), the left 
side of Eqn. |S| vanishes. Furthermore, when integrated over 
the entirety of an isolated system, S becomes zero and Eqn. 
|5] takes its popular form 



2T + W = G. 



(9) 



It is commonly taken for granted that W represents the 
gravitational binding energy of the system but it must be 
emphasised that this is strictly true only when computed 
over the totality of an isolated system. In general, W is not 
equivalent to the gravitational binding energy within r. 

In Figs. I15l and ll6l we present the temporal dep endence 
of the virial disequilibrium parameter iRickeil ll99^ for the 
gaseous and dark matter components of our evolving merger 
remnants following taccrete (when the system can be consid- 
ered to be a single remnant). (Although the baryonic com- 
ponent's contribution to the virial parameter of the total 
system is negligible, it is interesting to see if the observable 
gaseous component virialises on the same timescales as the 
dark matter component.) This virial parameter is defined as 
follows: 



2T 



W + S 



(10) 



where T, W and S are time averages computed to account 
for the fact that the canonical form of the virial theorem pre- 
sented in Eqn. |5] holds only when the quantities in Eqn. |S] 
are averaged over an interval in which the system is periodic. 
This is done using a dynamical time tdyn ~ ^/STr/lGGApc 
(where A is the over density being considered) for quantities 
involving the coUisionless component and the local isother- 
mal sound crossing time tsc = ^rUp/kT for quantities 
involving the gaseous component. Although the system may 
not be strictly periodic over these timescales, we have con- 
firmed that the term involving I in Eqn. El contributes neg- 
ligibly during the period we present in figures 1151 and 1161 
Following RSOl, we select the moment of virialisation to be 
when the system sustains a value of l^^j < 0.02. 

Overall, we find that the virial parameters for the gas 
and the dark matter track each other's general trends al- 
though the oscillations in the gas virial parameter tend to 
be of slightly lower amplitude and damp more quickly. As 
a result, the gas component typically relaxes (i.e. the gas 
virial parameter drops to and remains at |V| < 0.02) earlier 
than the dark matter component. There is not much differ- 
ence in the relaxation time of gas within R500 and -R200. On 
the other hand, the dark matter distribution within 7?5oo 
relaxes earlier than the full cluster. 

In more detail, our 1:1 mergers all virialise by treiax. 
The gaseous component generally virialises ^ 0.5 Gyr before 
the dark matter component which virialises almost exactly 
at treiax in all cases. Thus, the system is virialised approx- 
imately 4.5-5.5 Gyr (or 3-4 dynamical times) after the sec- 
ondary initially crosses -R200 of the primary (in good agree- 
ment the findings of previous studies, e.g. _Rocttieor ct all 
^23). The virialisation of our 3:1 mergers proceeds similarly 
to our 1:1 mergers for the dark matter. In the 10:1 case, 
the gaseous component within -R500 achieve virialisation at 
or well before our 1:1 and 3:1 mergers: immediately after 
taccrete iu the off-axis cases and well before in the head-on 
case. Within -R200, however, our 10:1 gaseous components 
reach virialisation much later (~ 2 Gyr) than our 1:1 and 
3:1 cases. The dark matter behaves similarly. Within iisoo 
the system is virialised at or before treiax while within i?200, 
the system never reaches our virialisation criteria. In these 
cases, the system sustains a state where Vd oscillates with 
an amplitude of ~ 5% until the end of our simulations. The 
retarded onset or absence of virialisation within ii200 for our 
10:1 mergers is a result of late infall of material dispersed 
beyond -R200 during the merger. 



5.2 Hydrostatic equilibrium 

Another independent measure of the equilibrium status of 
our merger remnants is the dynamical criterion of hydro- 
static equilibrium. For a pressure supported medium in a 
gravitational potential, hydrostatic equilibrium implies 



_idP _ GM{< r) 



dr 



(11) 
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Figure 15. Time evolution of the virial parameter computed over a volume of radius i?5oo (thick red) and i?200 (blue) for the gaseous 
component of our simulations. The evolution of the gas virial parameters for these two volumes is very similar. The horizontal dashed 
lines indicate our 2% virialisation criteria. The vertical dotted lines indicate taccrete a^nd t^^i^^ for each case. If we define the system 
to be relaxed when the virial parameter drops to and remains at \V\ < 0.02, then the gas component typically relaxes halfway between 

^accreteand ^reZax^* 



Given that X-ray derived mass determinations of clus- 
ters assume that this condition holds, its validity following 
mergers is of significant importance. Recent studies have 
found evidence suggesting that clusters for which hydro- 
static equilibrium can be accurately assumed are rare but 
no consideration has yet been given to whether systems can 
fail to be in hydrostatic equilibrium while appearing undis- 
turbed. In Fig. 1171 we plot the evolving hydrostatic disequi- 
librium parameter of our merger remnants following taccrete 
given by 



-1 dP 



^GM{< r) 



(12) 



The evolution of hydrostatic disequilibrium for our sys- 
tems all follow the same qualitative pattern. Shortly before 
taccrete, at both i?5oo and -R200, they experience significant 
positive excursions in H as the secondary core accretes to 
the remnant core. Following this, H is seen to oscillate with 
damped amplitude in each case with the oscillations at ii200 
lagging in phase behind those at 7?5oo by ~ 0.5 Gyr. The 
phase of these oscillations is generally such that at J?200, 
a period of positive H lasting ~ 2 Gyrwith an amplitude of 
10 — 20%, occurs at or shortly after taccrete- At R^o, shghtly 
lower amplitude (typically ~ 10 — 20%) and less coherent os- 
cillations persist until the end of our simulations. 

Regular disruptions of hydrostatic equilibrium such as 
those depicted in Fig. 1171 could translate into discrepancies 
in X-ray mass measurements obtained from methods which 



assume the va lidity of Eq n.[m We shall examine this issue in 
more detail in lPoole et a l. (2006b) where we will discuss the 
effects of mergers on global cluster observables and scaling 
relations. 



5.3 Comparison of apparent and formal states of 
relaxation 

In Section 2] we found that our merger remnants generally 
appear relaxed through centroid (our preferred method) and 
power ratio measures once the isophotes of our simulated 
2: = 0.1 Chandra observations appear undisturbed at treiax- 
In Table 121 we list the times at which our systems are 
formally virialised {t^iriai', measured for the dark matter) 
or in hydrostatic equilibrium {thydro', taken to be when the 
system sustains \H\ < 0.15) at i?5oo- Our systems are gen- 
erally virialised (both the gas and dark matter) within i?5oo 
at or near treiax while late accretion of material renders the 
system somewhat less virialised (5% versus our 2% viriali- 
sation criterion) at -R200 in our 10:1 mergers at treiax- We 
also list the redshifts at which the secondary would have to 
cross -R200 of the primary in order to be seen in each phase 
at redshifts of z — 0.1 and z = 0.5. We can see from this 
how unlikely virialised or hydrostatic clusters should be at 
z = 0.5. 

We thus find that our merger remnants are generally 
virialised at or shortly (< 2 Gyr) after appearing relaxed to 
visual inspection. However, even at this time merger rem- 
nants are adiabatically oscillating, exchanging energy be- 
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Figure 16. Same as Fig. ll5l but depicting the virialisation of the dark matter component of our simulations. The dark matter component 
typically relaxes slightly later than the gas component and the oscillations in V damp out more slowly. Also, the central regions of the 
systems relax faster than the full cluster^ . 
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Figure 17. Hydrostatic disequilibrium parameter measured at -R500 (thick red) and -R200 (blue). The horizontal dashed lines indicate 
15% equilibrium (the criteria used for thydro Table |^. The vertical dotted lines indicate t accrete and t^^i^x for each caset. 
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tween thermal and kinetic components. They are in hydro- 
static equilibrium only at the 10-20% level at this time and 
sustain this level of hydrostatic disequilibrium with little 
change after appearing relaxed. 



6 TRANSIENT STRUCTURES 

Our simulations reveal that the phenomena generally re- 
ferred to as "cold fronts" consist of several classes of tran- 
sient structure formed by the disruption of compact cool 
cores in off-axis mergers. Previous theoretical studies of cold- 
front production have focused primarily on the disruption of 
the secondary's core but our analysis reveals that both cores 
can be involved. Several other transient phenomena such as 
core surface brightness discontinuities and oscillations are 
naturally formed in the process as well. In this section we 
shall discuss in more detail the processes driving the for- 
mation of these structures, noting where possible the occur- 
rence of each in presently published observations of cluster 
mergers. 

6.1 Bridges 

At two points during a cluster merger, luminous bridges con- 
necting the interacting systems can be formed. In the first, 
the outer atmospheres of the merging clusters become com- 
pressed on the incident side of the collision during preinter- 
action. The resulting increase in gas density and tempera- 
ture between the systems leads to enhanced surface bright- 
ness between the two systems and the formation of a lumi- 
nous bridge. This effect is short-lived (~ 0.5 Gyr in the 1:1 
cases, less in others) but h as been observed in severa l sys- 
tems including A399/401 fSake lhou fc PonmanI liooi) and 
A1750 |Bc1so1c o^t al. 2005). 

As discussed in Section 13.31 material stripped by ram 
pressure from the interacting clusters becomes dispersed, 
forming comet-like tails which trail their surviving cores fol- 
lowing tciosest- Due to the density gradients of the systems 
they are passing through and the scaling of ram pressure 
with the density of the flowing material, this process is more 
efficient on the side of each core incident to the collision. As a 
result, overlapping enhancements of stripped material form 
behind each core, producing a luminous bridge which joins 
the two systems. 

Following tciosest, material from these regions moving 
in the plane of the orbit is driven together. It becomes par- 
ticularly concentrated in a curved planar feature oriented 
normal to the plane of the orbit. As a result, the surface 
brightness of this bridge is most bright and concentrated 
in ^-projections. In x-projections it is not discernible and 
in j/-projections it is significantly wider and more diffuse. 
The cores of both clusters appear disturbed at these times, 
distinguishing this process from the one acting during prein- 
teraction. 

The dependence of the brightness and shape of the 
bridge could be used as a useful cue when interpreting 
the orientations of obser ved merging systems. In A 1758 
jPavid fc Kempnerl |2004) for instance we see an excel- 
lent example of a system apparently exhibiting a post- 
interaction bridge. The two merging systems in this case 
(A1758N and A1758S) both have disturbed cores and 



there is an obvious bridge connecting them. The width 
of the bridge is more than 1.2 Mpc suggesting that the 
system may be significantly inclined to the plane of the 
sky. Ot her systems exhibitin g simil ar morphologies in clude 
A1644 iReiprich et alJl2004). A754 llHenrv et aljl2004) and 
A115 JGutierrez fc KrawczvnskilbOoiF ! 

6.2 Plumes 

In Section 13.31 we noted that in all of the off-axis cases we 
have studied, outwardly propagating collimated plumes of 
primary and secondary core material are produced between 
tciosest and tapo- Several mechanisms are involved in the 
production of these structures, a few of which have been 
discussed previously in the literature. iMathis et alJ (l2005h 
provide a good summary of theoretical research involving 
t he disruption of cold cores. Their work supports the claims 
of lHeinz et alJ (|20o3) who illustrate the important influence 
of ram pressure stripping during the disruption of merging 
cluster cores. However, they emphasise that other processes 
are likely essential to the formation of cold fronts. 

Our study confirms and expands upon this claim. Al- 
though ram pressure is certainly important for the disrup- 
tion of both clusters, it is primarily only responsible for the 
stripping of material roughly beyond i?2500 500 kpc (mod- 
erately high entropy material). It does not lead to the promi- 
nent low entropy radial plumes produced in our simulations. 
In our simulated Chandra images, the period during which 
ram pressure stripped material is bright enough to produce 
an observable "comet-like" surface brightness feature is brief 
(~ 0.5 Gyr) and occurs only at or shortly after tdosest- A dis- 
turbed merging system spends most of its time at or near 
apocentric passage. In our simulations, material stripped by 
ram pressure has dispersed and is no longer apparent in sur- 
face brightness by that point. Bright features of displaced 
core material are produced primarily by other mechanisms 
during this longer lived phase. 

In Fig. 1181 we illustrate the production of three low en- 
tropy plumes generated from the cores of the primary (Fig. 
118b .') and secondary fFig. ll8b 'l with a series of entropy maps 
covering the evolution of our 3:1 Vt/Vc ~ 0.15 simulation 
during first pericentric passage (Fig. 118b ) and apocentric 
passage (Fig. 118b '). In these figures, the secondary enters 
from the bottom left and passes underneath and near to the 
core of the primary. 

In Fig. 118b we see that two plume-like structures (la- 
belled "A" and "B") are produced from the primary core. 
The one labelled "A" is produced on the side roughly oppo- 
site from the point of closest contact with the secondary core 
at tciosest ■ This feature is generated from material which has 
become confined in a narrow structure by the convergence 
of the shock as it passes the primary core. This behaviour 
is produced by refraction of the shock due to the decline in 
the primary core's sound speed into the centre. 

The other two plumes ("B" and 'C' in Fig.lTSl are pro- 
duced through distinctly different means. In the off-axis sim- 
ulations we have studied, short-lived but strong negative ra- 
dial pressure gradients can be generated near the cores of the 
interacting clusters, leading to the ejection of low entropy 
material into the outer regions of the merger remnant. Since 
these gradients are often very localised, the ejection tends 
to be highly directional leading to the formation of large 
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Figure 18. The left panel of four plots shows entropy maps for our 3:1 vt/Vc = 0.15 simulation illustrating the ejection of two distinct 
"plumes" (labelled A and B) of low entropy material from the primary core in response to the pericentric passage of the secondary. 
The first is triggered by the shocks that precede the secondary, and the second by negative pressure gradients in the wake behind 
the secondary. The right panel of plots illustrates the ejection of a plume (labelled C) from the secondary at apocentre in the same 
simulation. The plumes eventually evolve into streams and most of the ejected low entropy material flows along the streams into the 
primary core. Each frame above is 1.5 Mpc on a side and is a z-projection of a 0.5 Mpc thick sliced. 



outwardly flowing coUimated plumes. Since the process is 
adiabatic, the ejected material retains its low entropy and 
eventually reaccretes to the core of the remnant. This oc- 
curs as low entropy streams of high velocity material over 
an extended period lasting well after taccrete (see Section 
lOl . 

Two mechanisms generally conspire to make this hap- 
pen: the sudden release of compressive forces due to changes 
in the relative radial motion of the systems and the pres- 
ence of low pressure regions in the direction toward which 
material is dispersed as a result. After pericentric passage, 
the primary's core is rapidly compressed and the wake of 
the secondary's motion creates a large low pressure region 
on the incident side of the collision. Low entropy material 
from the primary's core becomes ejected back along the sec- 
ondary's path as a result, forming plume "B"' which reaches 
a radius of ~ 800 kpc in this case. 

The secondary core is subject to these mechanisms at 
tapo when it turns around and returns to the primary core. 
The decline of ram pressure at that time results in the re- 
lease of compressive forces and the strong negative pressure 
gradient of the primary at such radii provides a low pressure 
region for the secondary core's material to disperse towards. 
These mechanisms are assisted by a gravitational "sling- 
shot" effect resulting from the lag of the sec ondary's gas 
core b ehind its dark matter halo, as noted bv lBialek et alJ 
i2002h . This lag is due to the additional pressure experi- 
enced by the gas but not the dark matter. The net result of 
these processes is the displacement of a significant portion 
of the surviving secondary core's material outwards to its 
tidal radius. This material is then efficiently stripped away 
by tidal effects forming plume 'C depicted in Fig. 118b . As a 
result, the plume formed from the secondary after apocen- 



tric passage assumes the morphology of a trailing tidal tail 
jToomre fc ToomrdlToT^) . More specifically, it has a shape 
which is curved in the same sense (from tail to core) as the 
orbit (see Fig. I18b 'l. As the system passes apocentre, this 
structure precesses until it is oriented radially with the pri- 
mary core. Plumes ejected from primary cores at tdosest are 
confined by the pressure gradients of the secondary's wake 
and do not share this shape nor evolution. 

The symmetry of our 1:1 off-axis mergers means that all 
of the processes discussed above act on both cores equally to 
create a pair of symmetric plumes. They connect at taccrete 
forming an "integral" shaped structure. Its curvature is 
much more pronounced in the vt/Vc ~ 0.4 case. Compar- 
ing Figs. |3to|S| we can see that the plumes formed from 
the secondary cores are quite similar between the 1:1 and 
3:1 cases. However, in the 3:1 cases the plumes formed from 
the primary core reach much smaller radii. 

The compressive forces generated in our 10:1 simula- 
tions are far less significant than in the 1:1 or 3:1 cases 
resulting in a much less significant plume of primary core 
material following tdosest- Furthermore, the resulting pres- 
sure gradients are far less localised and result in a broader 
and less coUimated displacement of material. By tapo, all 
of the material of the secondary beyond 300 kpc has been 
completely stripped and dispersed into the outer regions of 
the system. Part of the core survives until taccrete though, 
generating a small but undetectable plume after tapo- 

Several examples of observed systems exhibiting these 
morphologies exi st including A1758 JPavid fc Kemprieil 
2004) and A2744 iKempner fc Davidll20oir These systems 
possess disturbed accreting cores with curved shapes. In 
both cases, the authors suggest that they are products of 
ram pressure stripping but this interpretation is not consis- 
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tent with our simulations in either case: their large projected 
separations from their primary cores and their shape suggest 
they have ejected plumes similar to plume 'C. A1758N ap- 
pears to have passed in front of A1758S from east to west, 
transited apocentric passage and is ejecting a plume of ma- 
terial as it returns back toward A1758S. The A1758N plume 
is curved in the sense expected from such an orbit. A plume 
from the core of A1758S is apparent but lacks the same 
curved appearance. This may mean that A1758S is the more 
massive of the systems and has a plume such as plume "B" 
depicted in Fig. ll8h .. or that projection effects are distorting 
its appearance. 

A more secure example of a plume being ejected from 
a massi ve primary system su ch as plume "B" can be seen in 
A1644 iReiprich et alJ [20041) . As noted above, this system 
has a significant bridge joining the interacting systems. It 
appears to be at or near apocentric passage and the bridge's 
lack of curvature, the absence of curved plumes and the 
thickness of the bridge suggest that the system's orbit is 
highly inclined. Although the temperature map is cut-off 
near the main system to the south west, there is a strong 
suggestion of an elongated cold region similar to plume "B" 
seen in a y-projection. As we would expect from such an 
orientation, this feature is nearly coUinear with the bridge. 

Another good example of a system which appears to 
have had a plu me of primary core material ejected by a 
merger is A2744. lKempner fc Da^ j2004) suggest that this 
system is a 4:1 merger and should compare well to our 3:1 
off-axis cases. We conclude from the structure of the bridge 
joining the interacting systems and the shape of the sec- 
ondary's plume that this system is near apocentric passage 
and significantly inclined. The left-handed orientation of the 
secondary plume's curvature suggests that it has passed be- 
hind the primary in a north western direction and is now 
travelling towards us near apocentric passage. Consistent 
with this interpretation we find a cold region to the south 
east of the primary, where the secondary would have entered 
the system, ejecting a plume from the primary's core. 



6.3 Induced core rotation 

Due to the initial gas density gradients of the systems, the 
ram pressure acting on the cores as they pass each other 
is differential in off-axis cases. This induces torques on the 
gaseous components of each system causing them to rotate. 
Generally, the result is two co-rotating cores lasting from 
tciosest to the second pericentric passage at taccrete when 
the coherence of this motion is disturbed. Remarkably, the 
release of the pressure gradients which lead to plume "B" in 
Fig. ll8h applies enough torque on the primary core in our 3:1 
Vt/Vc = 0.15 case to cause the core to counter-rotate. This is 
the only instance in which we witness such counter rotation 
and more exploration of orbital and mass ratio parameter 
space is needed to understand how often this occurs. 

As a result of its induced rotation prior to apocentric 
passage, our secondary cores have prograde rotations when 
they eject the material which form the tidal tail plumes dis- 
cussed above. This enhances the efficiency of the process. 




Figure 19. Different time snapshots of the radial velocity maps 
(mass weighted) for the gaseous component in the 3:1 vt/Vc = 
0.15 simulation showing the collision of a high velocity stream (la- 
belled B) with the remnant core. In this specific case, the primary 
and the secondary cores merged at taccrete= 2.5 Gyr. Stream B 
consists of low entropy material that was originally part of the 
secondary core but shorn off during the course of the merger. 
Also shown is a stream feature (labelled A) that is not readily 
apparent in the surface brightness map but which corresponds 
to a stream associated with returning material originally ejected 
from the primary core during the creation of plume 'A' (see Fig. 
118b '). Each frame here is 3.0 Mpc on a side and is a z-projcction 
of a 0.5 Mpc thick sliced. 



6.4 Streams 

The plumes of material removed from the cores of merg- 
ing clusters discussed in Section 16.21 collect into large col- 
limated structures. Because of the reduced entropy of this 
material, its buoyancy is low and it ultimately accretes onto 
the remnant core after the remainder of the secondary core 
does so at taccrete- As a result, long lived (~ 2 Gyr) high 
velocity {v ~ 1000 kms~^) structures form in the cores of 
the remnant before the system relaxes. We distinguish these 
structures from plumes because they are inwardly propagat- 
ing high velocity structures of low (generally undetectable) 
surface brightness and are not associated with an apparent 
merging core (the secondary having been accreted and its 
observable presence destroyed at tacerete)- 

In Fig. 1191 we present a series of radial velocity maps 
for our 3:1 vt/Vc = 0.15 simulation illustrating the interac- 
tion of two low entropy streams (labelled "A" and "B" ) with 
the remnant core. The secondary core accretes from the top 
right at taccrete= 2.5 Gyr in this case and the stream formed 
from its disruption (labelled "B") follows it from the same 
direction. Another inwardly radial (but unapparent in sur- 
face brightness) feature (labelled "A") can be seen on the 
opposite side of the system in these frames. It is a long last- 
ing stream formed from plume "A" identified in Fig. 118b ,. 

Large scale bulk motions with velocities similar to these 
streams were observed in the cores of merger remnants in the 
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cosmological AMR simulations of iMotl et"all J2004ll . These 
authors raise concerns that such motions preclude the ex- 
istence of actively cooling cores. In our simulations, high 
velocity structures in the core remain highly coUimated and 
completely disrupt cooling of the core in only one case (3:1 
Vt/Vc = 0.15). Even in this extreme case, cooling resumes 
on rapid (~ 300 Myr) timescales once the stream finishes its 
accretion. This takes ~ 2Gyr. In all other 1:1 and 3:1 off- 
axis cases we have studied, the high velocity streams possess 
enough angular momentum when they accrete to miss the 
region of active cooling. In the 10:1 cases, the stream does 
not significantly disrupt cooling even when it does inter - 
act directly with the remnant core. In iPoole et all ll2006cll 
we will examine in more detail the effect these streams can 
have on the final structure of a merger remnant. 

The low entropy streams formed in the 1:1 off-axis cases 
are highly coUimated and accrete quickly soon after t accrete- 
In all cases, some portion of the streams are visible until 
tacerete whcu their appcarauce becomes disrupted by the 
second interaction of the cores. In the vt/Vc = 0.4 case, 
the accreting streams form a large low entropy spiral which 
evolves into a ring roughly 1 Mpc in diameter by treiax (see 
Fig-HJ- This structure then slowly dissolves, disappearing 
~ 2 Gyr after forming. 

The secondary cores in our 3:1 mergers become more 
significantly disrupted after tdosest than in the 1:1 cases. 
As a result, material stripped from the secondary core be- 
comes dispersed over a larger volume. Furthermore, it ac- 
cretes slightly later than tacerete and over a longer period. 
The streams are barely observable in our simulated Chandra 
observations for only a ~ 300 Myr following tacerete- 

Streams such as these have low surface brightness and 
are best observed in the temperature structure of clusters. 
They should manifest as cold fronts oriented radially to the 
primary core with no obvious secondary core being present. 
There are several examples of observed syst ems like this. In 
A576 for instance. iKempner fc DavidI ll2004f) note a "finger" 
of cool gas which extend s to the northwest of the system 
from its core. In A2255 iSakelliou fc PonmanllioO^ there 
is a distinct cool region to the southeast of the cluster core 
which is not associated with an obvious peak in X-ray bright- 
ness. Finally, in A2034 jKempner et al.l2003f) there is a large 
diffuse tail of relatively cool gas extending towards the south 
of the cluster centre with no obvious X-ray peak associated 
with it. 

6.5 Edges 

The enhanced spatial resolution of Chandra has revealed 
that the central regions of several clusters exhibit "edges" 
in their surface brightness as a result of large and localised 
gas density gradients. Among those arc RX J 1720. 1-1-2638 
(Mazzotta ct al. 2001), A496 (Dupkc fc W hite 2Qo3) and 
A 1795 llMarkevitch. Vikhhnin. fc Mazzottall200ll) . Systems 
with several of these features have a lso been found in- 
cluding A576 JKempner fc Davidll2004^ which has at least 
2 within 50kpc and ZW3146 (iForman et all 1200^ which 
has 3 at r ~ 20, 45 & 170 kpc. Attempts to account 
for these features have generally invoked bulk motions of 
cold core material with respect to its surrounding ICM. 
Explored mechanisms capable of driving suc h motions in- 
clude the activity of jets from central AGN d Fabian et alJ 
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Figure 20. Surface brightness maps of our 3:1 vt/Vc = 0.15 
simulation showing the formation of 2 "edges" associated with 
the accretion of the secondary (labelled A) and the stream of 
low entropy material that follows the secondary (labelled B). The 
first "edge" is weak feature that moves outward in radius and 
disappears soon after tacerete- The second "edge" is more well- 
defined and remains stationary during the accretion of the stream 
material, a process that can last for ~ 0.5 Gyr. We do not see any 
examples of "edges" produced by core oscillations, nor do we find 
evidence for significant core oscillations in our simulations. Each 
frame shown is a ^-projection 0.5 Mpc on a sidel. 



l2005fl or gravitati onally induced motions due to small pass- 
ing substructur e (iMark evitch. Vikhlinin. fc Mazzottall2boil : 
iTittlev fc HenriksenlbOOsF ! 

Occasionally when the streams discussed in Section 
accrete to the centre of the remnant, an interface forms be- 
tween it and the stationary core. Significant and abrupt 
jumps in density, temperature and pressure can result ac- 
companied by a jump in surface brightness reminiscent of 
observed edges. 

In Fig. 1201 we show the evolution of two such edges fol- 
lowing tacerete ioT OUT 3:1 Vt/Vc = 0.15 simulatiou. As the 
secondary core accretes to the primary core in this case, it 
does so on the top right side. This leads to an expanding 
edge (labelled "A") which starts at a radius of 120 kpc and 
disappears 0.3 Gyr later at a radius of 210 kpc. The stream 
which follows the accreting secondary core gets deflected to- 
wards the bottom of the core by this structure, forming a 
more signiflcant edge (labelled "B" ) at a radius of 40 kpc 
which lasts for 0.5 Gyr. 

In Section 16.41 we noted several systems which show 
evidence of streams of radially coUimated material. In two 
of these, A576 and A2034, there are obvious edge features 
present. In the case of A576, these edges are aligned prop- 
erly with the orientation of the stream to be a produced by 
it. In the case of A2034, the edges are apparently on the 
opposite side from the stream. In our simulations, multiple 
streams are typically produced by a single merger. It may 
be that an unseen stream (possibly produced by the same 
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event which generated the observed stream) is responsible 
for the observed edge in this case. It may also be that the 
observed stream is passing in front of or behind the core 
before making contact with it. 

No obvious edges are produced in our simulations by 
core oscillations. Edges produced through the accretion of 
streams are distinct, long lived (~ 500 Myr) and can yield 
multiple edges (such as those seen in A576) simultaneously. 
The last of these properties is difficult to account for through 
core oscillations but naturally explained as a product of ac- 
creting streams. 



7 SUMMARY AND CONCLUSIONS 

• This paper is the first in a series aimed at elucidat- 
ing the impact of mergers on relaxed X-ray clusters. For 
this study, we have carried out and analysed a suite of 
nine SPH simulations (incorporating the effects of cooling 
as well as star formation and associated SNe feedback) of 
two-body mergers between idealised clusters covering inter- 
esting and cosmologically relevant range of mass ratios and 
orbital properties. Our merging clusters are initialised with 
realistic initial dark matter and gas distributions, the latter 
resembling those seen in relaxed clusters with cool compact 
cores observed by Chandra and XMM. 

• Following a detailed qualitative analysis, we have iden- 
tified a generic progression for typical cluster mergers. The 
stages of this progression are: (i) pre-interaction, (ii) initial 
closest approach, (iii) apocentric passage, (iv) accretion and 
disruption of the secondary core and (v) relaxation. To fa- 
cilitate comparison and interpretation of observations, we 
describe in detail the interactions associated with the gas 
components of the two merging systems during each of the 
stages and quantify the times at which these stages arise 
relative to the beginning of the merger. We also show gas 
density, X-ray surface brightness, gas temperature, SZ and 
gas entropy maps at different points during the mergers. 

• Moderate and massive mergers tend to run for 
4.5-5.5 Gyr, while minor mergers last for up to ~ 
2 Gyr longer (in good agreement with the trend identified 
bv lTormen et al1l2004ft . This duration is the elapsed time 
from the point when the secondary's centre of mass crosses 
7?200 of the primary (the beginning) to the time when the 
gas distribution in the final system has relaxed and nearly 
returned to virial/hydrostatic equilibrium (the end). Typi- 
cally, off-axis mergers run longer than head-on mergers. 

• In simulated 50 ks Chandra observations (assuming z = 
0.1), all easily identifiable substructure generated by our 
merger events disappear during the secondary's second peri- 
centric passage of the primary core. However, the system 
takes an additional ~ 2 Gyr for its isophots to appear re- 
laxed. Temperature fiuctuations are not necessarily a reli- 
able indicator of the system's dynamical state. 

• We find that following the merger, the resultant system 
settles into virial equilibrium sooner than into hydrostatic 
equilibrium. In many cases, the final cluster exhibits oscil- 
lations that result in instantaneous deviations from hydro- 
static equilibrium by as much as 15% for as long as 7 Gyrs 
after the beginning of the merger. We also find that gener- 
ically complex patterns of entropy and temperature fiuctu- 
ations, at the level of AT/T ~ 20%, can persist in the final 



systems well past the point when the systems appears re- 
laxed. Based on our findings, we suggest that temperature 
fiuctuations of this order, without any other corroborating 
signatures of ongoing merger, should not be interpreted as 
evidence of a highly disturbed system. 

• We find that in none of the cases considered are the 
initial cool compact cores of the primary and the secondary 
destroyed during the course of the mergers. Instead, the two 
remnant cores eventually combine to form a new core that, 
depending on the final mass of the remnant, can have a 
greater cooling efficiency than either of its progenitors. We 
will discuss this behaviour in more detail in a forthcoming 
paper. 

• While the large scale coherent displacements between 
the dark matter and gas components in the merging system 
are not uncommon, they do not trigger any obvious surface 
brightness discontinuities. Within the central 50 kpc, the gas 
and the dark matter remain tightly coupled to each other. 
We do not see any evidence of core oscillations or of any 
associated heating. 

• We investigated the efficacy of three measures used to 
quantify the degree to which observed clusters are disturbed. 
These are (i) centroid variance, (ii) power ratios, and (iii) 
X-ray surface brightness/projected mass displacement. We 
find that the centroid variance is the best of the three. It best 
refiects the state of the cluster as determined by visual ex- 
amination. It also provides an excellent indicator of how far 
the system is from virial and hydrostatic equilibrium. The 
power ratios are very sensitive to noise and once this is fac- 
tored in, they provide an unambiguous signal only when the 
system is highly disturbed. The X-ray/mass displacement 
are extremely sensitive to weak but long-lived perturbations 
in the dark matter distribution caused by low mass mergers. 

• During the course of the merger, we find a variety of 
transient features arising and disappearing in the projected 
temperature, entropy and surface brightness maps, most of 
which resemble "cold fronts" observed in recent Chandra 
and XMM observations. While "cold fronts" are generally 
associated with the disruption of the secondary's core, we 
find that they can be caused by several different mech- 
anisms. The resulting features can potentially be distin- 
guished based on their morphological properties. Interested 
readers are referred to Section [H] for examples of clusters 
whose recently published X-ray images show examples of 
the different types of transient structures. We propose the 
following classification scheme based on their origin and ap- 
pearance: 

- "Comet-like" tails: Ram pressure strips the outer at- 
mospheres of secondary systems forming very brief (~ 
0.5 Gyr) comet-like morphologies for both cluster cores 
(except the primary in 10:1 mergers). This morphology is 
observable only between first pericentric and apocentric 
passages. 

- Bridges: Ram pressure disperses core material from 
both cluster cores in 1:1 and 3:1 cases. A luminous 
moderate-entropy bridge resulting from the convergence 
of the two resulting regions of dispersed material is pro- 
duced. This structure remains observable in 50 ks Chandra 
exposures at a = 0.1 until secondary core accretion in our 
1:1 mergers and until apocentric passage in our 3:1 merg- 
ers. The morphology of this structure could be useful as 
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a cue for determining the inclination of observed merger 
systems. 

- Plumes: Pressure gradients produced by the sec- 
ondary's motion at tdosest and tapo can lead to the ejec- 
tion of significant amounts of the primary and secondary 
cores' gas. This material forms large coUimated plumes 
of low-entropy material. At apocentric passage, material 
ejected in this way from the secondary's core becomes 
tidally stripped and adopts the morphology of a trail- 
ing tidal tail. The majority of observed disturbed merging 
cores likely owe their appearance to this process, rather 
than ram pressure stripping. 

- Streams: Low entropy streams form from these 
plumes and accrete to the centre of the remnant following 
the second pericentric passage and disruption of the sec- 
ondary core at t accrete- This material remains coUimated 
and typically has enough angular momentum to avoid a 
direct impa ct with the remna nt core. This mitigates the 
concerns of iMotl et alJ (|20o3) who suggest that high ve- 
locity structures formed in mergers may prevent the for- 
mation of compact cool cores. These streams may also 
play an important role in fuelling the activity of central 
AGN in clusters. 

- Induced core rotation: Off-axis mergers subject the 
gaseous cores of both merging clusters to torques which 
induce rotation. In one simulation (our 3:1 Vt/Ve ~ 0.15 
case), the resulting rotation of the primary's core is retro- 
grade to the secondary's orbit. In all cases, the coherence 
of this rotation is completely disrupted by the disruption 
of the cores during second pericentric passage. 

- Edges: Streams can accrete directly onto the core of 
the remnant, generating strong "edges" in surface bright- 
ness and gradients in temperature and density. The result- 
ing features are similar to several cold front systems whose 
morphology had previously been attributed to AGN ac- 
tivity or motion of the gaseous core within its dark matter 
potential. 

We would like to reiterate that our simulations do not 
incorporate the physical processes associated with magnetic 
fields, pressure from cosmic rays or thermal conduction nor 
do we adequately model the effects of turbulence. The influ- 
ence of these processes on clusters remain contentious issues 
and we can not make any confident predictions at this point 
regarding their effects on our results. 
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